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Verf. überträgt seine T’heorie der numerischen Kubatur !), die inzwischen in anderen Arbeiten ?) ?) abgewan- 
delt worden ist, auf den Fall allgemeiner Integrationsbereiche im k-dimensionalen Raum und gibt einige 
Anwendungsbeispiele. 


In 1936, the author gave a theory of numerical cubature*) that has, since then, been differently treated in 
other papers ?)®). The iheory is generalized to an arbitrary range of integration in the k-dimensional space, 
and some examples are given. 


L’auleur applique sa theorie de cubature numerique!) ‚qui, depuis ce temps, a &i& modifiee dans d’auires 
iravaux2)®), au cas de spheres d’integration generales dans l’espace de la dimension k et donne quelques. 
exemples pratiques. 


ABTOp mepeHocHT CBOW TeOPHIO WHCAIeHHOU Kydarypst 1936 r.!), BIOCHIEeNCTBUNH BUNOHS- 
MeHeHHyIO B APyrux paborax?)®), Ha cıyyaä OÖMMHX KOHTPYpoB WHTETPAAIOB B k-MepHOoM 
upocTpaucrtse. JlaHo HecKOAbEO IIPHMepOoB IIPHMeHCHHN. 


Der entscheidende Gedanke meiner Theorie besteht darin, daß die Integration nicht 
auf einer vorhergehenden Interpolation der zu integrierenden Funktion aufgebaut wird. In 
der Tat ist die Aufgabe, eine Funktion in brauchbarer Weise zu interpolieren, viel weiter 
reichend und viel schwieriger als die, ein bestimmtes Integral der Funktion zu berechnen. Diese 
Loslösung des Quadraturproblems von dem der Interpolation (im eindimensionalen Fall) 
tritt, soweit ich sehe, zum erstenmal in einer Arbeit von W. Wirtinger ?) zutage, die in 
dieser Zeitschrift 1933 erschienen ist. Eine eingehende Studie des eindimensionalen Falles, die 
auch Stieltjes- Integrale und mehrfach gegebene Punkte umfaßt, habe ich 1935 an anderer 
Stelle) veröffentlicht. 

Die folgenden Rechnungen. sind vollkommen elementar. Sie erfordern kein anderes 
Hilfsmittel als das der Produktintegration in einer Dimension. Die Einfachheit des Resul- 
tats wird evident, wenn man einen Blick auf die Hauptformel (Il) in Abschn. 5 wirft. Diese 
Gleichung kann nicht verbessert werden. Dagegen bleibt manches zu tun, um in besonderen 
Fällen bessere, d.h. speziellen Bedingungen angepaßte, Abschätzungen für die Fehlergrenze 
zu finden. 


1. Voraussetzungen und Bezeichnungen 


Es sei eine stetige, beschränkte Funktion f von k Variablen ®%,, &, .... %, in einem Gebiet 
@ des k-dimensionalen Raumes gegeben, über das sie zu integrieren ist. Von fsetzen wir voraus, 
daß die Funktion in jedem Punkte von @ nach jeder Richtung mindestens einmal, im all- 
gemeinen m-mal stetig differenzierbar ist. Von dem beschränkten Gebiet G nehmen wir folgen- 
des an. Es gibt einen Punkt O im Innern von @ von der Beschaffenheit, daß jeder Halbstrahl 
durch O die Kontur von @ genau einmal trifft (Sternbereich). Ist @ konvex, so tut jeder innere 
Punkt den Dienst. Ist ein Gebiet von einer endlichen Anzahl analytischer Flächen begrenzt, 
so läßt es sich immer in endlich viele Sternbereiche zerlegen. Unsere Formeln lassen sich also 
in allen solchen Fällen verwenden. Man könnte aber auch, wenn man einige Komplikation der 
Rechnung in Kauf nimmt, die Theorie direkt, d.h. ohne Zerlegung, für andere als Stern- 
bereiche brauchbar machen. 

Die Werte der Variablen in O, die wir auch die Koordinaten von O nennen, seien %, = En 
%,= &,...% = &. Ein Halbstrahl durch O, also eine Richtung im k-dimensionalen Raum, 
wird durch (k — 1) Verhältniszahlen (#, —.&) : (2 — &) :...(8 —&,) oder durch (k —1) 
geeignet gewählte Winkelgrößen @, d, y,... bestimmt. Ist r der Abstand eines beliebigen 
Punktes von O, so kann die Kontur des Gebietes @, daher auch dieses selbst durch eine Glei- 


chung der Form 
Bel, W,e.«) Q Er (2) 


1) Revue mathömatique de l’Union Interbalkanique, t. 1, 1936, fasc. 1. 

2) MauroPicone, Annali della Scuola Normale Superiore di Pisa. Serie III. Vol. V. Fasc. III—IV, 1951. 
®) Carlo Birindelli, Compositio Mathematica. Vol. 10, 1952, S. 117—167. 

“4 W Wirtinger, Z.angew. Math. Mech. 13, 1933, S. 166—168. 

5) R. v. Mises, J.reine u. angew. Mathem. 174, 1935, S. 56—67. 
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Er 


olumelement des Kegels, ı 
die Größe rt-1dr dS. Z.B. im =" 
Radius 1 ist, wird dS = dg, nämlich das Differe 
ment“ ist r dr dp. Für k = 3, bei der üblichen Ww; 
dS= sin ddp dd und das Volumelement wird r 
Kenntnis dieser Ausdrücke wird im weiteren nicht benöti 
in dem zu berechnenden Integral % 
Lane 


wo dV ursprünglich für das Produkt da, da,...da; steht (und 
solcher Zeichen vertritt) anstelle von dV auch i 


1 dr dS 
gesetzt werden kann. Die Definition (2) nimmt so die Form an: 


= [ffr dr is = | [rear] ds. 


Hierbei geschieht die Integration über dr bei konstant gehaltener Richtung und die darauf E 
folgende Integration über d$ ist auf die gesamte Einheitskugel auszudehnen. 


2. Problemstellung 

Wir nehmen an, daß in n Punkten P,, P,,... P„, die im Innern oder am Rande von @ 
liegen®) (und von denen zunächst keiner mit O zusammenfallen soll) die Werte der Funktion f 
bekannt seien. Wir nennen den Funktionswert im v-ten Punkte f,; der Abstand dieses Punktes 
von O heiße r, und die Koordinaten, d.h. die Werte der unabhängigen Variablen in P, seien 
= 123,3...Kk;v=1,2,...n). Jedem Punkt wird eine ‚‚Gewichtszahl‘“ A, zugeordnet ° 
und die Summe der Produkte A,f,wirdalsein Näherungswert des Integrals J 
betrachtet. Die Differenz 


DI -SAh= [FI ZAf, SE 


bezeichnet man als den Fehler der Integrationsformel. Unsere Aufgabe ist, D zu berechnen 
und aus dem Ausdruck für D Gesichtspunkte für die Wah! der Punkte P, und der Gewichte A, 
abzuleiten. Wenn wir für Deinen Ausdruck finden, der nur von den Ableitungen m-ter Ordnung 
von f abhängt und der verschwindet, sobald alle diese Ableitungen in allen Punkten von @ null 
sind, so sagen wir, wir hätten eine Integrationsformel m-ter Ordnung vor 
uns. 

Eine Formel m-ter Ordnung muß jedenfalls die Bedingung erfüllen, daß D —=0, sobald 
für fein beliebiges Polynom von niederem als m-tem Grad gesetzt wird. Alle solche Polynome 
sind enthalten unter den linearen Kombinationen der Potenzprodukte 


&. &ü & 
u a 


wenn &1, &, . . . & alle natürlichen Zahlen, einschließlich der Null, sind mit einer Summe kleiner 
als m. Die Zahl der von einander verschiedenen Produkte dieser Art ist 


fm +k—1\ (m+k-—1)! 
Zu —| n En ee) 


Die Wahl der Punkte P,. P,,... P„ und der Gewichte A, 4, ... A, wird also für eine 
Formel m-ter Ordnung jedenfalls durch Z,,, „ Bedingungsgleichungen (man. wird später sehen, 
daß diese auch hinreichend sind) eingeschränkt. Werden die Punkte willkürlich gewählt, so 
muß man die Gewichte so bestimmen, daß sie Z,,z Gleichungen erfüllen. Es sind also im all- 
gemeinen mindestens Z,,; Punkte für eine Formel der Ordnung m erforderlich. Es kann aber 
vorkommen, daß bei besonderer Wahl der n Punkte einzelne der A, gleich Null werden, so daß 
tatsächlich weniger Punkte in Gl. (4) auftreten. Eine Formel m-ter Ordnung, in der weniger als 
die durch (5) bestimmte Anzahl Zm,» von Punkten verwendet wird, nennen wir eine „‚„Ver- 


°) Es ist hier angenommen, daß die n Punkte alle von einander verschieden sind. D 
Punkte ist in der Fußn. 5) zitierten Arbeit behandelt. Man 
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feinerte“ Formel m-ter Ordnung. Beispiele verfeinerter Formeln im eindimensionalen Raum 
sind die Simpsonsche Formel mitn—=3, m —4(während Z,,, —4) und die Formeln der 
Gaußschen Quadraturmethode, wo n bei geradem m gleich m/2, sonst gleich (m +1)/2 ist 
bei Fa = m. Alle diese Fälle werden von unserer Theorie mit umfaßt, ebenso wie solche mit 
N > Zum, 

Zur Einführung in die allgemeine Aufgabe stellen wir eine Besprechung der Formel erster 
Ordnung voran. F 


3. Integrationsformel erster Ordnung 
Von der Funktion f setzen wir jetzt nur voraus, daß sie stetige Ableitungen erster 


Ordnung in jeder Richtung, d. h. nach jeder Variablen &,, %,, . . . &;, besitzt. In einem beliebigen 
Punkt (außer O) heiße f’ die Ableitung von fin der Richtung OP, also 
er U — of 
I ne 0 ar 0%, = pr 9300, #5 re (6) 


Ist P,einer der gewählten Punkte und Pein Punkt auf dem Halbstrahl OP, in der Entfernung 
r von O, so nennen wir f, (r) den Wert der Ableitung in Richtung OP, im Punkte P. Es ist also 
J» (r,) die radiale Ableitung in P,, wofür wir auch kürzer f/ schreiben, und es ist f/ (0) die Ab- 
leitung in O in der Richtung OP,. Andere Ableitungen in O werden für die Formel erster 
Ordnung nicht benutzt. Es ist aber für das Spätere wichtig, folgendes zu bemerken. 

Wir bezeichen mit f, den Wert von fin O und mit /) die Ableitung von fin O für eine 
beliebige Richtung OP. Ist r der Abstand des Punktes P von O und sind &,, %, . . . %, seine 
Koordinaten, so gilt 

; 0 ) 0 ? 
rfo (5) + 4 + (de 2 ee) 
wobei für 9f/0x,; die Ableitung nach x; im Punkte O einzusetzen ist. Nach (6’) ist rf, eine 
lineare Funktion der Koordinaten &,, &, .. . %. 

Von den Gewichtszahlen setzen wir für die Formel erster Ordnung nur voraus, daß ihre 

Summe gleich dem Volumen V des Integrationsbereiches ist: 


Arts AA eV a a Ze el) 
Das Volumen kann zufolge von Gl. (3) in folgenden Formen geschrieben werden: 


R 
var [| [mals [ mas. N 5: 
ü 0 


Nunmehr beweisen wir den Hauptsatz’), der später auf den Fall der Formeln 
höherer Ordnung ausgedehnt werden wird: 

Ist fveinmalistetig differenzierbarund)genügen dierA,der 
Bedingung (7), so gilt: 


r 


Di av — 2 Afı— [| moral-34 ne, ee... 


M)=L(R-r), NMian)=1. 


Um dies zu beweisen, wenden wir auf das erste Integral über dr auf der rechten Seite 
Produktintegration an. Es ist wegen M,(R)=0 und Mi(r)= — r!=! 


R R 
| M,(r)f’ = —tRep+ | ifar. 
Ö ö 


Da ferner 
1, 
HM: N,(r, t,)f»(r) dr =f, — 
0 
?) Die Ableitungen in diesem Abschnitt sind umständlicher als für den unmittelbaren Zweck nötig wäre; 


sie sind so gefaßt, wie es für die spätere Verallgemeinerung wünschenswert ist. 
14* 
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so wird die rechte Seite von (I) | 
1 : n n { 
Pen, j RS + | | \ fm dr) B-EAfthEA. (9) 


iedi j j ite Glied 
Das erste Glied ist zufolge (8) gleich — Vf, das letzte zufolge (7) gleich 14 Jo. Das zwei 
ist gemäß (3) gleich J, demnach der ganze Ausdruck (9) identisch mit der linken Seite von (I) 
w.z. b.w. ’ 
Um (I) für eine Integrationsformel auszunutzen, beachten wir, daß M;(r) im ganzen 
Intervall (0, R) nicht negativ ist und daß 


R R 
| mo =, = u (10). 

Demnach hat das erste Glied auf der rechten Seite von (I) die obere "Grenze 
Fax [|] 7% ar] N HAN A Be ee. ür: 

und man kann aus (I) schließen 

IDI< Kine | SraP + rl] ee. il; 
Sind die A, nicht negativ, was in der Regel der Fall ist, so kann man mehr aussagen, nämlich 
San 7 U —faas Zr Ar Da fine fm Zrrhr >. > (12) 


In der Gestalt (11) oder (12) wird später das allgemeine Resultat für Formeln der m-ten 
Ordnung ausgesprochen werden. Hier, bei m — 1, können wir noch einen andern Schluß ziehen, 
indem wir für das letzte Glied in (I) mit Rücksicht auf N, (r, r,)—=1 den gleichwertigen Aus- 
druck 3 A,(f, —f,) setzen. Dann wird wegen (7) 


jsav-pr=s| Sms arlas 
und daher nach (10) und (10°) 


a} 1 Fe 


Dies gibt eine Abschätzung für den Unterschied eines über V genommenen Mittelwertes von f 
und dem Wert von fin O, wobei für O jeder Punkt genommen werden darf, von dem aus jeder 
Halbstrahl die Kontur des Gebietes nur einmal trifft, also jeder innere Punkt im Falle eines 
konvexen Gebietes. Die beste Annäherung bietet im allgemeinen jener Punkt, für den das 
Polarmoment ein Minimum ist. Im Falle einer Kugel um O mit dem Radius R (im dreidimen- 
sionalen Raum) ergibt (13) 


3 E il 3 F 
ZRfmasg [FAV-SS ZT Rfne EN 


Andere Beispiele werden später (Abschn. 7) erörtert werden. 


4. Vorbereitung für den allgemeinen Hauptsatz 

a) An Stelle der ersten radialen Ableitung f’ von f werden wir jetzt die Ableitung u-ter 
Ordnung f“ benutzen, das ist, fin einem beliebigen Punkt P (außer O) u-mal nach r bei fest- 
gehaltener Richtung OP differenziert. Den Wert von f® in einem Punkte P, der in einem Ab- 
stand r von O auf dem Halbstrahl OP, liegt, bezeichnen wir mit F (r) und schreiben für f (r,) 
kürzer f{9. Der Wert von f(# (r) für r — 0, also (0) ist die w-te Ableitung inO in der Rich- 
tung OP.,. 

Die radiale Ableitung f(* läßt sich in bekannter Weise durch die u-ten Ableitungen nach 


% “, ©, ausdrücken. In der üblichen Schreibweise des symbolischen Produktes ist, analog 
zu 
MI % 7 50 % — 50 aa 
J Fe Sl MO er 70% £ el, 


Re 
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Die w-te Ableitung von f im Punkte O in einer beliebigen Richtung OP, wobei P die Koordi- 
naten %,, %,, ***, 2,und den Abstand r von O besitzt, ist analog zu (6’) gegeben durch 


pa - ++ nr (15) 
ON von 2 20%, k Ko RT D 


wenn hier die partiellen Ableitungen von fnach den &,, 2,, *** z,im Punkte O eingesetzt wer- 
den. Es ist also r#f ein Polynom u-ten Grades in den Koordinaten ©, 2%, *** 2. 

b) An Stelle der Größen M,, N,, die in (I) eingeführt wurden, werden wir Polynome M,, 
N, verwenden, die folgendermaßen definiert sind 


M,,(r) -[M (o)doe, M,;,(r) (Rt erh) ee FR): 
N,.(r, r,) = [N (rd N lerne, SA (17). 


Es ist ohne weiteres zu sehen, daß M,(r) ein Polynom vom Grade (k +u — 1) ist, das für 
r —=R verschwindet und sonst im Intervall von 0 bis R nur positive Werte annimmt. Ebenso 
ist N,(r, r,) vom Grade (u— 1), verschwindet (außer bei u =1) an der Stelle r =r, und ist 
sonst positiv im Intervall von 0 bis r,. 

Für die N bestätigt man unmittelbar die Formel 


1 
a Tl Ne N ee LEE 


Für die M gilt die Darstellung 
R 


are 
Muh (r) = M J Imre 
0 


yk+u 
FEN ar: 


Denn dies stimmt für u= 0, und Differentiation der rechten Seite nach (— r) stellt den gleichen 
Ausdruck wieder her, nur mit u ersetzt durch (u — 1). 

Für später brauchen wir noch die Werte der M und N an der Steller =0. Für N gibt 
(18) und (17) | 


"u 
Ve j Heyne a; a (20), 
während für M aus (19) folgt ; 
1.ReTae 
er 
c) An Stelle der einen Gl. (7) werden jetzt die Gewichte A,, 4,, ..- A, dem ganzen System 
von Gleichungen unterworfen, die besagen, daß die Formel 


SpAd,=jfaV 


korrektist, wenn feine Konstante oder ein Polynom in %,, %, --- 2, vom ersten, zweiten, 
dritten, ... Grad bezeichnet. Es seien %1,, &2,, -- + %%, die Koordinaten (Variablenwerte) des 
Punktes P,®=1,2,... n). Dann haben wir 


R R 
M „1 (0) = j M,,(e) de= N NEN 
() 


N 
1 Gleichung _ 0-ter Ordnung ISA, y 


v1 


k Gleichungen 1-ter EEE ENTE N ae 
v1 
n ( (22). 
” 2-ter „ I AeaeV; BA... k 
‚1 


ie 


ns 3-ter £ 


v 


0,2, = [0,0 0V;, u A=l,2,...% 
usw. ) 


Die Anzahl aller Gleichungen bis zur Ordnung (m — 1) einschließlich ist in (5) angegeben wor- 
den. Die Koeffizienten der A, in (22) hängen nur von der Wahl der n Punkte P, ab, die rechten 
Seiten nur von der Gestalt des Integrationsgebietes. Man sieht leicht ein, daß die Gleichungen 
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Bedingungen ausdrücken, die von der Wahl des Koordinaten-Systems unabhängig sind, d.h. 
unverändert bleiben, wenn die x, durch irgend welche lineare Funktionen der x, ersetzt werden. 
Z.B. besagen die Gleichungen nullter und erster Ordnung zusammen, daß Massen von den 
Größen A,, Ay --- An, diein den Punkten P, konzentriert sind, denselben Massenmittelpunkt 
besitzen wie eine über das ganze Gebiet @ gleichförmig verteilte Masse von der Größe ihrer 
Summe. Nimmt man die Gleichungen zweiter Ordnung hinzu, so bedeutet dies, daß die beiden 
Massensysteme auch in den Trägheitsmomenten übereinstimmen usw. 


Wir werden die Gl. (22) als de Momentenbedingungen bezeichnen. 


5. Der allgemeine Hauptsatz 


Nach den voranstehenden Vorbereitungen wird es leicht sein, folgendes zu beweisen: 

Hauptsatz. Wenn die Funktion fim ganzen Integrationsge- 
biet stetige Ableitungenin jeder Richtung bis zur Ordnungm 
besitzt und die Gewichtszahlen A, 4A, --- 4, den Momentenbe- 
dingungen (22) bis einschließlich zur Ordnung (m— I) genügen, 
Solgıl tabume 1,2, 3a 


D= Jrav - 8 a5= [| Morwa]as— ZAafn.an) 1a... m), 


wobei die MundN die durch (16) und (17) definierten Funktionen 
sind. 

Die Behauptung ist für m =1 in Abschn. 3 bewiesen worden. Es ist also nur noch zu 
zeigen, daß wenn sie für irgend ein u zurecht besteht, sie auch für («+-1) gilt, vorausgesetzt, 
daß (u-+-1) nicht größer als mist. Mit anderen Worten, es ist zu beweisen, daß die rechte Seite 
von (II), falls u < m, ihren Wert nicht ändert, wenn darin u durch (u +1) ersetzt wird. 

Nennen wir zur Abkürzung K, und Z, die beiden Integrale über dr in (II), so daß die 
rechte Seite lautet 

N 
[K,dS — 3 AL, 
v1 


so ist der zu beweisende Satz 


[Kun K)S-3 Aland)... (2). 
Nun hat man 
R * 
Ru Ma fd Made 2... (28). 


Da aber nach der Definition (16) M, die negative erste Ableitung von M,+1, ist der Ausdruck 
in der Klammer die Ableitung des Produktes M,+1 f® und da alle M für r— R verschwinden, 
ergibt (25) 

Kun — Ku = — Mr). 
Hier setzen wir für den ersten Faktor rechts den letzten der in (21) ermittelten Ausdrücke und 


erhalten 


R 

Ku E— K, == fin e ykt+u—1 dy 
Ba 0 

und schließlich für die linke Seite von (24) 


1 
[am mas || mtwaras =, [mnwar ..... 0 


zufolge der in (3) erklärten Gleichheiten. 
In ähnlicher Weise verfahren wir mit der rechten Seite von (24). Es ist 


r, 


Lu+ == Hr =) Nana dr, 


0 


und da wi i i : : \ r 
da wieder N, die negativ genommene Ableitung von N,+ı ist und dieses für r — r, Ver- 
schwindet, so wird 


Lay — u = — Ny41(0, 7.) [@ (0,7,) = — Im fi (0, 5) 
ul 7 


a  ı 
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wobei der Wert für N,+1(0, r,) fürr —=0 aus (18) genommen ist. Die rechte Seite von (24) hat 
mithin den Wert 


n n 
2 Allan) =— a SA er 


Unsere Aufgabe ist also darauf zurückgeführt zu zeigen, daß die Ausdrücke (26) und (27) 
gleichen Wert haben, d.i. daß 


n 
SAT) move er 
v=1 


Im vorhergehenden Abschnitt ist unter a)gezeigt worden, daß r#f{® ein Polynom u-ten Grades 
in %4, %9, ««. @,ist. Der Wert dieses Polynoms in einem Punkt P ist r* mal der u-ten Ableitung 
von fin O für die Richtung OP. Der Wert in P, ist gleich r“ mal der Ableitung in O für 
Richtung OP,, also r} mal f(#(0, r,). Somit drückt (28) nur die Tatsache aus, daß für ein ge- 
wisses Polynom u-ten Grades das Integral über dV gleich ist der Summe der mit den Gewichts- 
zahlen A, multiplizierten Polynomwerte in den Punkten P,. Dau <m — 1, ist dieses Faktum 
durch die Forderung, daß die A, den Momentenbedingungen bis zur Ordnung (m — 1) genügen, 
sichergestellt. Damit ist (28) und rückläufig (24) und schließlich der Hauptsatz (II) selbst in 
vollem Umfang bewiesen. 


6. Die allgemeine Integrationsformeln 
Nachdem wir für Din (II) einen genauen Ausdruck gefunden haben, ist es nicht schwer, 
Abschätzungen für Dzu geben. Es kommt uns dabei zustatten, daß die Funktionen M und N, 
wie in Abschn. 4 gezeigt, in ihren Integrationsintervallen 0, R bzw. 0, r nur positive Werte an- 
nehmen. Es folgt also jedenfalls aus (II) 


Diet | [|| 90 ir) ds +2 1A,| In r,) i Be ui), 


Das hier auftretende Integral von M,(r) dr ist in (21) berechnet worden. Nehmen wir den 
letzten der Ausdrücke (21), so wird das Doppelintegral in (29) gleich 


R 
= [rar a3. [mar. 
ul |, u! 


Ebenso liefert die Definition von N in (18) und (16) 
T, R 1 
[ IN lt r,) dr = Nur1(0, 2, = u! iz (31). 
4 | 


Demnach nimmt (29) die Form an 
DI< Ilm | [mrar+Flale| a Wet: 
° v—1l 


Die Berechnung einer oberen Schranke für |D| verlangt somit nur die Kenntnis des u-ten po- 
laren Moments des Gebietes @. Ist u eine gerade Zahl, so ist r* ein Polynom vom Grade win 
den z,, nämlich 

ln — Et ee... (83). 
Bei ungeradem u wird man in den meisten Fällen besser die Sch warzsche Ungleichheit be- 


nutzen, die besagt, daß 5 
SredV <yYVfr*dV a ER N 


Die Abschätzung (32) gilt für jedes ganzzahlige u bis einschließlich u — m. 

Zu Beginn von Abschn. 2 war angenommen worden, daß keiner der n Punkte P, mit dem 
Ursprung O zusammenfällt. Dies war für den Rechnungsgang wünschenswert. Da aber die 
Endformel (32) für beliebige Werte der Abstände r, gilt, muß sie auch richtig bleiben, wenn 
eines der r, gegen Null geht. Man darf also ohne weiteres auch zulassen, daß O einer der Punkte 
P, ist; der betreffende A,-Wert tritt dann in (32) nicht explizite auf. 

Ist nicht nur für den absoluten Betrag von f® eine Schranke gegeben, sondern eine obere 
und untere Schranke für f® selbst und sind die A, alle positiv, so kann man (32) verschärfen zu 


) 

D<ı rm [rar-1m Su e| 
u! MAT BT u 

1 ei El a 

Dell [rar 2 ar) 
u! mın mar 2 
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Mit min = — [mas geht (35) in (32) über. Man verbessert die Abschätzung, wenn man als Fak- 
toren der Integrale in (35) die Extremwerte nur derradialen (statt jeder) Ableitung u-ter 
Ordnung und als Faktoren der Summen den größten, bzw. kleinsten der Ableitungswerte auf 
den n Halbstrahlen nach den Punkten P, einsetzt. Doch wird die dazu nötige Detailkenntnis 
von fim allgemeinen schwer zu erlangen sein. 

Sind alle A, nicht-negativ, was in der Regel zutrifft, so bleibt als hauptsächlicher Mangel 
der Abschätzung (32) der, daß die Summe von zwei Gliedern in ihr erscheint, während die 
wahre Größenordnung des Fehlers ihre Differenzist. Im eindimensionalen Fall wird bei 
der Simpson schen Regel dieser Nachteil durch eine andere Art der Rechnung vermieden ?). 
Ob etwas derartiges auch bei der mehrdimensionalen Integration möglich ist, muß dahinge- 
stellt bleiben. Es ist zu bedenken, daß n Punkte auf einer Geraden (einem einzigen Halb- 
strahl) etwas anderes bedeuten als n Punkte in einem auch nur zweidimensionalen Gebiet. In 
den folgenden Beispielen werden wir uns nur der Abschätzungen (32) bzw. (35) bedienen. 


7. Einige Beispiele 
a)Kugelim k-dimensionalen Raum. Esseieine Funktion f von %k Varia- 
blen %,, &, -»- %, zuintegrieren über das Gebiet, in dem 
HR Hti MER „ 2. she een 
Wir wählen einen der Punkte P, im Zentrum und geben ihm das Gewicht A,, ferner je zwei 
Punkte auf jeder der k-Achsen in Abständen + r, mit den Gewichten A,. Dann sind der Sym- 
metrie wegen mit den drei Gleichungen 
Asch = 2A = 424), 2kAın— [rtdV EN 
alle Momentenbedingungen bis einschließlich der Ordnung 5 befriedigt. Die Formeln, zu denen 
wir in dieser Weise gelangen, sind somit in hohem Maße ‚‚verfeinerte‘‘ Formeln im Sinne von 
Abschn. 2. Ist S die Oberfläche der Einheitskugel — Kenntnis der Größe von 8 ist für die 
Rechnung nicht erforderlich — so gilt nach (3) - 


R R 
7-8 [ma -s% [rar=s | ra, R2V |rar= tz Rev 
3 BEN k +2 i k +4 
() \ (38). 
Daher liefern die Gl. (37) 
k+2 k +4 4 
ee Eee Br} — 2m: 
rl Tarp bh ein De 35) % 
und die Fehlergrenze nach (32) wird 
14 @% k(k +2) 
rl et 
DI<|f Imaz 6; E+6t ra a N 


Für den Fall k=2, Kreisgebiet in der Ebene, 5 Punkte, findet sich mit V = R?x aus (39) 
und (40) 


= oR, Ai R?r, A=Z R’r, 
DI < mas gg ir 
Für k = 3, Kugel im dreidimensionalen Raum, 7 Punkte wird mit V = - Rz 
A= IR, A=:Rn, A= Rn, 
[DIS Inne a Po. 


Natürlich kann man bei derselben Wahl der Punkte und Gewichte auch Fehlerformeln niederer 
Ordnung benutzen, wenn die Abschätzung der (radialen) Ableitung sechster Ordnung zu schwie- 
rig erscheint. 

b) Einpunktformel zweiter Ordnung. In Abschn. 3 ist der Fall, daß 
nur ein Funktionswert für die Auswertung des Integrals benutzt wird, behandelt worden. 


®) Daß auch die übliche Feblerformel der Simpsonschen Regel nicht di i 

‚’) Daß iblic] gel nicht die genaue Fehlergrenze lief. 
habe ich in einer Note in dieser Zeitschrift gezeigt: Z. angew. Math. Mech. 13, 1933, S. 53—56. De Re 
handelte Problem ist für eine Verallgemeinerung auf mehrdimensionale Räume zu schwierig. 
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Die Wahl des Punktes im Innern des Integrationsgebietes war im wesentlichen willkürlich und 
die Fehlerformel war von erster Ordnung. Eine in den meisten Beispielen günstigere Lösung 
erhält man, wenn man den Schwerpunkt von @als den Punkt wählt, dessen f-Wert ver- 
. wendet wird. Es sind dann die Momentenbedingungen bis zur ersten Ordnung erfüllt und man 
kann (32) mit u =2 anwenden. 

Nehmen wir an, daß der Schwerpunkt auch zugleich als Pol O dienen kann, so folgt aus (32) 


D< |" Ins 4 | rrav le a I IE ER LEN 


Diese Formel kann in der Weise Verwendung finden, daß das gegebene Integrationsgebiet zu- 
nächst in eine Anzahl von Teilen 


Ven+n+t--V, 


zerlegt wird. Sind dann fj, fa --- fs die Funktionswerte in den s Einzelschwerpunkten und 
01> 0% - ++ Or die polaren Trägheitsradien dieser Teilgebiete, so hat man 
1 0 | fi a Br vum 
W753 I Y7 = If Imaz 372% n= If Imaz 2 Oma *» * * =* (42). 


Eine Zerlegung in rechtwinklige Zellen (und eventuell Teilzellen am Rande) bietet sich als ein- 
faches Beispiel dar. Aber auch Dreiecke (Simplex)-Netze können Verwendung finden. 

c) Rechtwinklige Zelle Simpson-Formel. Ein Rechteck (oder Par- 
allelogramm) läßt sich leicht in ein Quadrat transformieren; das Entsprechende gilt von den 
analogen Gebilden im k-dimensionalen Raum. Es genügt also den Fall des ‚, Würfels“ in k Di- 
mensionen zu betrachten. 

In Analogie zu der Simpsonschen Formel in einer Dimension bilden wir eine (ver- 
feinerte) Formel vierter Ordnung, indem wir dem Mittelpunkt des Würfels ein Gewicht A, und 
jedem Eckpunkt ein Gewicht A, zuschreiben. Aus der k-fachen Symmetrie folgt dann, daß alle 
Momentenbedingungen ungerader Ordnung erfüllt sind. Ist « die Seitenlänge des Würfels, so 
lauten die Bedingungen nullter und zweiter Ordnung 


2 k+2 
At A, = at, #48 | war, 
also 
2 | 
A=za, Amt za. 


Bezeichnen wir mit f, den Wert von fim Mittelpunkt des Würfels und mit 2 das arithmetische 
Mittel der 2#*-Werte in den Ecken, so wird 


1 = 
A zetah +.) 3 


Die Fehlergrenze rechnet sich mit Hilfe von 


k k(k—1 
[rar=r| Ay ta 1) [adv = ar] + = 
zu k( k) 
1 (1+5% 
IV — ak+4 
IDIS If" Imaz z @ 190 °’ 
so daß a. 
1 = 2a trdı Kill 0%) 
ZI R TI / q4 |fIV Ku en > 
Y / 3 | mas = 


Für k = 1 ist die durch (43) gegebene Fehlergrenze 4mal größer als die für die eindimensionale 
Simpsonsche Regel auf anderm Wege abgeleitete. 
Wird die Seitenlänge des Würfels in s gleiche Teile geteilt und die voranstehende Über- 
legung auf jeden Teilwürfel angewandt, so wird das Fehlerglied rechts in (43) s*-mal kleiner. 
In der Integrationsformel auf der linken Seite von (43) erscheint jeder der s* f-Werte, der dem 


2 EHRSER 
Mittelpunkt eines der kleinen Würfel entspricht mit dem Faktor E3 s-=. Die übrigen (s+1)* 


Werte zerfallen in (k +1) Kategorien, von denen die x-te (x = 1,2, ... k+1) aus 


k 4140-191 
„ie 11H, 
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Einzelwerten besteht, deren jeder das Gewicht 
il 


= Yan—1—k s* 


erhält. Die Summe der Gewichtszahlen in (43) wird dann 
k+1 


2 4 1 k lc — 
2 ER ee 
wie es sein muß. 


Mehrere andere Beispiele für den Fall k=2 sind in der in Fußnote 1 zitierten Arbeit zu 
finden. 
Eingegangen am 26. Januar 1953. 


Über Iterationsverfahren für Gleichungen und 
Gleichungssysteme. (I. Teil) 
Von Rudolf Ludwig in Braunschweig !) 


Unier dem Gesichtspunkt der Ordnung eines Iterationsverfahrens ( Bodewig-Hartree) gewinnt man durch 
einen allgemeinen Ansatz mit willkürlichen Funktionen eine große Anzahl von I terationsverfahren, unter denen 
sich bekannte wie das Newtonsche befinden; aber auch neue und solche, die noch willkürliche Funktionen zur 
Konvergenzverbesserung enthalten. 


The paper deals with iteration methods and considers them from the point of view of ihe order of the 
process ( Bodewig-Hartree). A great number of methods is obtained by a statement containing undelermined 
functions. Among ihem, there is the well-known Newton’s method as well as some new that contain arbitrary 
functions which may be used to improve Ihe convergence. 


Sous le point de vue de l’ordre d’un procede d’iteration ( Bodewig-Hariree) on obtient un grand nombre 
de methodes d’iteration par un arrangement general de fonctions indelerminees. Parmi ces möthodes se 
trouvent telles connues que celle de Newton, mais aussi de nouvelles et telles quelles contiennent encore des 
fonctions arbitraires applicables au perfectionnement de la convergence. 


B cTaTbe paccMaTPuBamTcA UTEPANHOHHLIe METONEI C TOYKH 3peHHA TOPANKA IIPoMecca 
nrepannn (Bonesura-Xaprpn). IlocpencTBoMm HeKOTOpoH OÖmmeii HCXONHOH CXeM&EI C IIPOH3BOJIB- 
HEIMH DyHEIMAMH HOAY4YAETCH ÖOMBIIO2 KOAHYECTBO HTEPALHOHHEIX METONOB, K YHCHY KOTOPBIX 
HPUHANIEIKaT H USBECTHEBIe, KAK HAIIPHMeEPp, cIIOCOd HBOTOHA, a TAKHE H HOBEI®, H TAKHE, KOTOPEI® 
eme cONeP:KaT IIPOM3BONIBHLIE DYHKIIUH MIA YIyYIIeHHA CXONHMOCTH. 


Einleitung 

Es sollen hier JV 2) für Gleichungen und in einem später folgenden Teil II für Gleichungs- 
systeme behandelt werden. Die Betrachtungen gehen aus von dem ‚„Konvergenzverhalten‘‘, 
auf das wohl zuerst Bode wig [2]°) ausführlich hingewiesen hat, und auf den von Har- 
tree [4] geprägten Begriff der ‚‚Ordnung eines JV‘““. In sehr allgemeiner Weise werden JV 
verschiedener Ordnung durch einen Ansatz mit willkürlichen Funktionen hergeleitet. Aus der 
Vielfalt der verschiedenen möglichen Verfahren kristallisieren sich bekannte Verfahren wie das 
Newto.nsche als einfache Fälle heraus. Aber auch andere Verfahren, die noch willkürliche 
Funktionen enthalten, sind wichtig und vorteilhaft. Wie die abschließenden Betrachtungen 
zeigen, nimmt der Rechenaufwand für einen Schritt mit der Ordnung zu, aber es sinkt für die 
geforderte Genauigkeit die Zahl der notwendigen Schritte. Der praktische Rechner kann dar- 
a: je nach der vorliegenden Funktion und ihren Ableitungen, sich den günstigsten Fall aus- 
wählen. 

1. Iterationsverfahren für Gleichungen mit einer Unbekannten 


11.Ein Konvergenzmaßfür Iterationsverfahren 
Es soll ausgegangen werden von einem JV in der Form 


J bedeutet hier eine im allgemeinen reelle Funktion einer Variablen x, aber fund x können auch 
Vektoren bedeuten (s. Teil II). Davon, daß man fauch als Operator auffassen kann, soll hier 
kein Gebrauch gemacht werden. Damit (1,1) als JV einen Sinn hat, muß die Konvergenz nach- 
gewiesen werden. Für den praktischen Rechner genügt die nachgewiesene Konvergenz allein 


!) Diese Arbeit wurde von der Mathematisch-Naturwissenschaftlichen Fakultät der Techniseben Hoch- 
schule Braunschweig als Habilitationsschrift genehmigt. Der Verf. ist den heiden Berichtern, Herrn Prof. Dr. F. 
Rehboc kund Herrn Prof. Dr. L. Collatz für wertvolle Ratschläge sehr zu Dank verpflichtet. 

2) JV = Iterationsverfahren; JV/]=JV1l. Ordnung, JVı = JV 2. Ordnung, JVk = JV k-ter Ordnung. 

8) 5. Literaturangaben am Ende der Arbeit, 
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nicht, denn für numerische Rechnungen darf die Konvergenz nicht so langsam sein, daß der 
Rechenaufwand unangemessen groß wird. 
Zur Erklärung der Konvergenzverhältnisse definiert Bode wig[2]: 
Eine konvergente Folge &,, &, **, &, mit dem Grenzwert heißtkonvergent 
vom Grade %, wenn 
A 
lim - 
Nn> © (%, => z)* 
dabei soll keine positive Zahl sein. Bezeichnet man die Abweichung der einzelnen Glieder der 
Folge vom Grenzwert mit 


le 0 fs le Wa te (1.2), 


a ER TRE Tre EN 
so kann man statt (1,2) auch schreiben: 
Iime-erl.cnn.n. et, 


N 
Entwickelt man r„}ı nach Potenzen von r, und beginnt diese Entwicklung mit er, so konver- 
giert die Folge vom Grade *. 
In gleicher Weise hat Hartree [4] de Ordnung eines JV eingeführt. Für die 
Gleichung 
Diey=f(2) te 0A en 
auf die Form (1,1) gebracht, führt man nach (1,3) die Abweichung des n-ten Iterationsschrittes 


von einer Lösung zein. Die Entwicklung von r„+1in eine Taylorreihe nach Potenzen von 
r„ ergibt: 


= 1 Pe 
ART ORH nei) 


Es wird dann definiert: 


Ein JV für eine Gleichung &,41 = f(x,) heißt von k-ter Ordnnng, 
wenn für eine Lösung & gilt: (1,6) 


Fo=f@)=:..= fm @)=0, aber, FOR)Z0. 


Diese Definition ist eine Spezialisierung der allgemeinen Definition von Bode wig für kon- 
vergente Folgen auf den Fall der JV für Gleichungen. 


E27 Spezielle Bedingungen für JPY, 
Für eine algebraische oder transzendente Gleichung 
D(x) = 0 
ist das Newto.nsche Verfahren weithin bekannt: 
D (2) 
In+1 — In D(x ) s 
n 
Die in (1,1) gegebene Form wird analog der Gestalt des Newtonschen Verfahrens in den 
folgenden Betrachtungen sehr oft die Gestalt annehmen: 


nr =) = In B(a,) ee ee: EEE. 


Zur Aufstellung der Bedingungen, die nach der Definition (1,6) die Funktionen A(&) und B(x) 
erfüllen müssen, damit (2,1) ein JV, darstellt, muß f (x,) an der Stelle nach Taylor ent- 
wickelt werden. Um an späterer Stelle die Ausführungen nicht zu unterbrechen, sollen schon 
jetzt diese Bedingungsgleichungen angegeben werden. Unter der Voraussetzung 


B@) #0. - EIN (2,2) 
lauten die Bedingungen 
aytürren Jr ll) Aß)=0; 
b) JVı: (I) A@)=0, A’) = B(e); 
c) JVy : IM) A@)=0, 4@=B@), 4'(@) = 2B(f); 
d) JVYıv: (IV) A@)=0, 4'@) = B(x), A’(&)=2B(@), A’) = 3 B’R); 
e) DNB ee, 1er, 1). 


4) Von den hier und in der Folge auftretenden Funktionen wird gefordert (sofern nicht ausdrücklich etwas 
anderes gesagt wird), daßsie und ibre Ableitungen bis zu einer genügend hohen Ordnung in einem abgeschlossenen 
Bereich reell und reellwertig, eindeutig, beschränkt und stetig sein sollen. 
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1,3. Iterationsverfahren zweiter bis vierter Ordnung 

Nach einer Methode, die Do m b [3] in ähnlicher Weise für algebraische Gleichungen an- 
gewandt hat, werden für Gleichungen jeder Art — algebraische und transzendente — durch 
Einführung willkürlicher Funktionen JV der Form @„41 = f(@,) möglichst großer Allgemein- 
heit hergeleitet. 

Für eine Gleichung 


B(x) =— 0 A a ee ter a u He (3,1) 
wird zunächst vorausgesetzt, daß sie sich auf die Form >= 
®(2)=:H() -G(e)=D Paaren (3,2) 


bringen läßt. Durch Einführung der beiden willkürlichen Funktionen »(x) und g(z) gibt man 


(3,2) die Gestalt 
« [p(z) H(&) + g(#)] = p(a) @(@) + wg(a). 
Diese Gleichung erhält die Gestalt (1,1) eines JV, indem man durch die Klammer der linken 
Seite dividiert: | 
nn Pla) &a) + ade) 2 POLE Le ee 
p(=) H(«) + q(@) p(«) H(z) + g(@) 


= r(a) dla)  _„_AR) 
Or ernee dr  2m a er 


Die rechte Seite des JV hat die allgemeine Form eines solchen: 
In = + nF (En) 


mit der passend gewählten Zahlenfolge c, 61°", m» "- 
Wegen der speziellen Gestalt (2,1) muß nach (2,2) für eine Lösung & vorausgesetzt werden: 


Be Dale) eo 20. ee (3,4). 
Die Bedingung (I) A(&) = 0 ist wegen ®(x) = 0 erfüllt, also stellt (3,3) für alle Funk- 
tionen p(&) und g(x) ein JV, dar (Konvergenz vorausgesetzt). 
Für ein JV,; muß weiterhin die zweite Bedingungsgleichung (II) erfüllt sein: 


4) — B@) = PR) DR) — p@) H@) — a@)=0, 


also 


also 


1%) z z 
„6 = PO) — H@ (3,5). 
Diese Bedingung schränkt die Funktionen p und gein und zwar dadurch, daß an der noch zu 
bestimmenden Stelle für das Verhältnis beider Funktionswerte ein bestimmter Zahlwert vor- 
geschrieben wird. 

Ahnliche Bedingungsgleichungen wie (3,5) werden in der Folge sehr häufig auftreten. 
Die Mannigfaltigkeit der Funktionen, die nach (3,5) zugelassen werden können, ist natürlich 
sehr groß, da (3,5) nur eine Aussage macht über den Funktionswert an einer Stelle & (oder evtl. 
auch über das Verhalten an allen Lösungsstellen), die allerdings nicht explizit bekannt ist, 
sondern als Lösung der Gleichung erst bestimmt werden soll. Daher wird man die zu bestim- 
menden Funktionen aus den gegebenen Funktionen der zu lösenden Gleichungen aufbauen. 
Durch die Einführung solcher willkürlicher Funktionen läßt sich das Konvergenzverhalten 
günstig beeinflussen, und es lassen sich viele in der Literatur zerstreut vorkommenden JV 
unter einem gemeinsamen Gesichtspunkt behandeln. 

Zur Bestimmung von p(z) und g(x) sollen folgende Möglichkeiten aufgezeigt werden: 

a) (3,5) werde identisch erfüllt für alle : 


ge) = Pla) (Be) Hin). 2 2.0 er 
In (3,3) eingesetzt folgt daraus die rechte Seite für das bekannte Ne wtonsche Verfahren: 


D(%) 


Ko=3— gel Van san le. 


Aus der Voraussetzung (3,4) folgt, daß nur dann ein JV,, vorliegt, wenn ®’(@) #0, d.h. & 
muß eine einfache Nullstelle sein. Bei einer mehrfachen Nullstelle liegt nur ein JV, vor. Dar- 
auf weist übrigens auch Bode wig [2]hin. 
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b) Geht man aus von der Bedingungsgleichung (3,5) in der Form 
4’(#) — B(a)=(p®/ —pyH—g=0, 
so folgt daraus ein JV mit der rechten Seite: 


x) D(x 
| fa)=»— ne ap RD ee (3,8), 


unter der Voraussetzung p(&) #0. Durch die willkürliche Funktion p(x) läßt sich entweder 
Tr Konvergenz für einen Anfangswert &, erreichen oder die Konvergenzgeschwindigkeit er- 
öhen. 


ce) Wählt man für die in (3,8) vorkommende freie Funktion p(e): 
p(x) = et”, k = const, 
so spezialisiert sich die rechte Seite des JV zu 


E D(e) 
U Eiler zer 


PER REES 


Eine vorteilhafte Wahl von % siehe (3,13). 
d) Weitere willkürliche Funktionen kann man dadurch einführen, daß man als Lösung der 
Funktionalgleichung (3,5) ansetzt: 


ne Four +Aı(a) Bla) + Asa) Ba) +... 


Durch Einsetzen in (3,3) kommt man schließlich auf die folgende sehr allgemeine Form eines 
Ve: 


Da) + Ale) Bela) HA) Dad) +: 
D’(x) + ua) D(x) + ul) D(a) + - 


mit den willkürlichen Funktionen 4; (2), u. (2), @= 1,2, °*:). 

e) Die willkürlichen Funktionen A;(#) und u;(&) können dazu dienen, Verfahren von 
mindestens 3. Ordnung festzulegen. Setzt man in (3,10) A,(#) =},(a)=-::-=0 
und u;(2) = u;(2) =: =(), geht also aus von 


-_ $(a) + Aa)B*(a) | 
Ta) De) EEE" 


so folgt aus der dritten Gleichung (III), indem diese Gleichung identisch erfüllt sein soll: 
2A) Da) 2) =EDlM) 2.22.2222... 


fs) =2— (Par = mel9sL0) 


Daraus läßt sich folgern: 
) u(e)=0, Aa) =PD’(x)/2 D°(x), also: 


Ss we ante 
pP) Aa) =0, ul) = — ®’(2)/2 D’?(x), also: 
2.0(2) © 

19-:- Te en Tee 


Diese beiden Verfahren dritter Ordnung sind sehr ähnlich; welches Verfahren rechentechnisch 
günstiger anzuwenden ist, hängt von der Funktion Dab (s. Beispiele). 
Durch Vergleich von (3,9) und (3,12) folgt, daß bei der Wahl von 
&’ (*) 
20 (=) 
5) Nach dieser Formel rechnet Do m b [3] für algebraische Gleichungen, ohne sie explizit anzugeben und 
ohne anscheinend den engen Zusammenhang zum Ne wt onschen Verfahren zu sehen. 


8) Dies ist das J Vı11, das meist im Zusammenhang mit dem Ne wtoon schen Verfahren aus der Ta ylor- 
reihe hergeleitet wird. 


k= 
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(3,9) von 3. Ordnung werden müßte. Da x noch unbekannt ist, kann man &% durch einen Nähe- 
rungswert x, ersetzen. Man erhält ein JV,,, das einem JVırr „‚benachbart‘‘ ist, da in der 
Taylorentwicklung (1,5) das erste Glied c,r» einen um so kleineren Koeffizienten hat, je 
besser der Näherungswert a,ist. Das JV lautet: 


x 2D(x) D(%,) 4. 
a 2 D’ (&,)®’ (&,) — D’ (2%) D(,) 


y) Aa)=k De) ‚ u(e) = (% = 5) ve k = const, also 


Men: 


02 100" © | 
an GV 


o: + 5) oo ? | 


Fa) = e— 


Ein Vergleich von (3,14) mit dem Ne wtonschen Verfahren als JV;, zeigt, daß dieses durch 
eine Korrekturfunktion der Form 

D(2) 

P(e) 

zu einem JV,7, verbessert wird. Es sei noch darauf hingewiesen, daß man mit einem Nähe- 
rungswert «,, den man konstant läßt, ein JV7r 


f(@) = — K(a) 


D(&) 
D' (2) 
erhalten kann, dessen Konvergenzverhalten günstiger sein kann. 

6) Für ein JVıv sind durch Gl. (*) die ersten drei Bedingungen (IV) bereits erfüllt. Aus 
der vierten folgt: x 

62(a) ®’ (=) 8" (2) — 3u(x) ©" (x) + 612 (2) D*(e) + w(z) B’(a)] = 28" (2). 

Nimmt man die nach « differenzierte Gl. (*) hinzu, so lassen sich A und w’ eliminieren und A 
und u bestimmen. Die rechte Seite des JV heißt dann: 


En+ı = En — K(%) 


f(a) er e 6 &'2 &'’ .- 2 &' &'' rn S &®'’2 ® 
= 768282..20'97.-0020.0, 


(IV 


Auch hier tritt eine solche Korrekturfunktion auf, durch die aus dem JV „ein JVıy wird. 


f) Aus den hergeleiteten JV sieht man, daß man die rechte Seite eines JV als Polynom in 

®(x) mit Koeffizientenfunktionen A; (2) ansetzen kann: 
Fa) = # + Aı(@) Pla) +Rz(a) D’(a) +A,(@) Dia) + 
Die Anwendung der Bedingungsgleichungen (D—(IV) ergibt: 
; “ 1 

Al) = 0, PP 5 06 lese 
Erfüllt man diese erste Bedingung näherungsweise, indem man & durch x, ersetzt, so folgt als 
rechte Seite des JV 


Ö 
So=3-gn (IV 
=(,i ı D(%) ’ 
An =0, = 3,4, A, ze 00, = lee 
2) Ü }1(8) De) A,(®) IX ZIEH const führt zu 
n Da) 
fe) =® ae Dia) \, (IV) 


Dieses Verfahren entspricht (3,13), ist aber, da }, nur näherungsweise festgelegt ist, von 
2. Ordnung. 


?) Daß diese oft als Vereinfachung des Ne wt on schen Verfahrens angegebene Form von 1. Ordnung ist, 
wird meist unzureichend begründet (z. B. [8]). 
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)A@)=0, 1-45; 


ee ne ie] 
D'(x) 2 D°%(x) P' (x) 


6D%(z) 
(3,15) verwandtes JVıy mit 


DD 
(2) = 030: 60% u ‚ (JVıv)®) ee (3,18). 


Zur Verringerung des Rechenaufwandes kann man A,(&,) = const wählen, erhält allerdings da- 
durch nur ein Jr: 


= 0) ergibt ein 


D(2) D’(a) De) Dia) [,D%a) un 
ano, Seren 


Io, (dV m) (3,19) 


Beispiel: Iterative Berechnung der n-ten Wurzel aus der Gleichung 
D(e)=ır —a—=0. 
Fürn=3,a=30, x, = 3ergibt sich: 
a) Nach Newton: 
Jo)=x 
b) Nach (3,11): 


ar — a .—ıa 


non’ an 9 111 111% (JV7)- 


a 1 +2 


c) Nach (3,15): 
ee el er ont un 2, Sm a 
2n+2) 0" +2(n— 1)a yarıı? 1 0 423 +5a 322 > ya! 
%, = 3107 235 


Der exakte Wurzelwert ist: x = 3,107 2325. 
Zusammenfassung: 
Für eine gegebene Gleichung Ö(xz) = 0 werden JV von der Form m+1ı = f(#.) 


aufgestellt. f(x) kann in vielfältiger Weise aus der Funktion ®(z), ihren Ablei- 
tungen und evtl. aus willkürlichen Funktionen p(&), A(z) usw. aufgebaut werden: 


J(a) = A[Pla)]. 
Die hergeleiteten Verfahren (3,7)—(3,19) mit Angabe der Ordnung und des erfor- 
derlichen Rechenaufwandes sind zusammengestellt in Tabelle I, S. 223£. 


(3,20) 


1,4. Konvergenzbedingungen 

Zum Vergleich der Verfahren verschiedener Ordnung werden solche Konvergenzbe- 
dingungen angegeben, die diesen Vergleich bequem zulassen. Für spezielle Verfahren, wie z. B. 
das Ne wto.n sche, lassen sich sehr oft schärfere Abschätzungen formulieren. 

Für ein JV, 


In-+1 —t(&) es Geste een ee a (4,1) 


mit einer Lösung, also &=f(&), folgt durch Ta ylorentwicklung nach (1,5): 
f®(E 
In+1 — Zr ie . . . . ee ee le un . (4,2). 


Die £, liegen in einem Bereich ®, der die Näherungswerte &, und die Wurzel x enthält. Die 
obere Schranke der k-ten Ableitung in diesem Bereich sei: 


ES M,. ol. „n—1; 


dann kann man abschätzen: 
kr_ı 


else A RE 


8) Dies JVıv wird meist aus der Ta ylor reihe als Erweiterung des Ne wt on sschen Verfahrens gefunden, 


A . . angew. Math. Mech. 
216 Ludwig, Iterationsverfahren für Gleichungen und Gleichungssysteme. I. Teil A en 
“ld. ı Ludwig, Nberationsveranren TUE TerSaEBs 


Daraus folgt die hinreichende Konvergenzbedingung: |[r,®"1M, <1, ausführlicher 
geschrieben: 


| Io, — | f® (E) man < kt| ERIRERANEN, ME 


Dabei ist &in B so zu wählen daß der Betrag der k-ten Ableitung seinen maximalen Wert 


» 


annimmt. Speziell folgen daraus die Konvergenzbedingungen für 
a) ein JV;: IF (E)Imas <1, 
b) ein Vz: | —Zllf"(Elmae< 21, 
c) ein Vr: 100 — ER] (E)lmas < 31. 

Es gilt also: 


Ein JV, konvergiert, wenn es einen Bereich ®* gibt, der x ent- 
hält und für den |f’(«)|<1. Für ein JV,(k 22) gibt es immer 
einen Anfangswert x, so, daß |a, — #1] f®(E) mas <k!, d.h. daß 
das Verfahren für diesen Anfangswert konvergiert. 


(4,5) 


1,5. Verbesserung einer Iterationsfolge 

Von einer Iterationsfolge werden die letzten drei aufeinander folgenden Werte betrachtet, 
die hier mit %,, %1, &, bezeichnet werden sollen. Aus diesen drei Werten soll nun ein verbesserter 
Wert & berechnet werden, ohne zu iterieren, also ohne Berechnung neuer Funktionswerte und 
deren Ableitungen. Für die Folge eines JV, erhält man dadurch eine Formel, die wohl zuerst 
von Steffensen[6]angegeben wurde, aber sie wird auch von Hartree[4]und Samuel- 
son [5] erwähnt. 

Nach der Definition eines JV, ist nach (1,5) 


1 a 
Ge OL Er ze rEan oT suntee 


Wird r, als ‚„‚genügend‘‘ klein vorausgesetzt, so kann man den ersten Summanden der rechten 
Seite als Näherungsausdruck betrachten: 


ß 1 = 
Na+ı ”4r rk mit (4 = rl (2) . 


Setzt man in dieser Näherung das Gleichheitszeichen, so muß der wahre Wert & durch den 
verbesserten Näherungswert & und r, durch 7,, a, durch a, ersetzt werden, also: 


n=&+ To; =itn=it+ N, n-itn=ärtat + (2). 
Die Elimination von a, und r, ergibt: 
(H— 2 Ele ar nee 
Unter Einführung des Differenzenschemas der Iterationsfolge: 
% = n,—% 
0 1 0 
3... 07, | Pe ee 
%; 
folgt aus (5,3) 
a 1 a 
(Kö, — 61) (a, — "+ 5 K(k +1) art 4... Hör =0 .. (5,4). 
Dies ist eine Gleichung k-ten Grades für &. Speziell folgt daraus 
a) für ein JV;: 9) m — 2) +F—=0, 
ee a 
= AAN Arne 


b) für ein JV: 
2-5) —- HIER - HM ÄI=-0 v0 . (5,6); 

c) für ein IV gr: 
(3 9 — 1) (to — 8)? + 6 dg (2, — ©)? + -M)HA=0.... (5,7). 


®) (8,5) wird auch von Willers [7]angegeben und bemerkt, wie sich spät i 
gangs-J V nicht zu konvergieren braucht. De a I Tr ee 


zn 
- 
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Schreibt man allgemein für ein JV, 
2=%, +&%6, : und N OR N ae (5,8), 

so erhält man die Korrekturen eg, aus den Gleichungen 

a) fürein J,: A Ya—y?=0; 

b) für ein JV/r: Q—yJ)&+(l Paar =; 

e) fürein Im: B- N) +3l1+2y—-WE+l+3Yy+3Pr 39), — Y=0. 
Das Bild 1 zeigt e, bis e,in Abhängigkeit von y aufgetragen. 

Bei der unter b)aufgeführten Korrektur & 


für ein J V;r ergibt die Ausrechnung EBRZEEBESFTDG 
3 ee] DRENSERZENER 
£ 1+4 
3-4 |, u. 69 


ER Yes i 
Dabei muß stets y> — — sein. Eine genauere 


Untersuchung lehrt, daß das Vorzeichen der 
Wurzel stets positiv sein muß. Bei der ku- 
bischen Gleichung für ein JV7ı, kann die Wur- 
zeldiskussion unterbleiben, da man praktisch 
diese Gleichung ebenfalls nur näherungsweise 
lösen wird und ein Näherungswert für & leicht 
anzugeben ist. 


BNZENERBEN-ENE 
LER ee 
NABNEIBBEETEIN 
L\. IR are 


Um einen Überblick über die Größe des Fehlers von & zu erhalten, kann man folgender- 
maßen vorgehen: 


Dem Gleichungssystem (5,1) der Näherungswerte soll das der exakten Werte gegenüber- 
gestellt werden. Es ist: 


n=:+n 
tt ale) (5,10) 
B—=L tr +41 (E,) str far 
= +". Hrn ar (Ed + ar (Er) + rar Et” 
Setzt man noch 
Nomen er, Aue ad N. (5,18) 


und zieht die entsprechenden Gleichungen der Systeme (5,1) und (5,11) voneinander ab, so 
erhält man: 


X+R,=0, 
X +rRad + LEER +... HD) RS, tra li), 
K tr Ayla 0) HI HT RR, att 
Hört [Eat tan) A (Sr rad (Ed +. Arkrtat,, 
+ Dale) artnet. 
Ersetzt man nach der Abspaltung der Differenzen 
04? a, ie NE, Ga nca re, 
so erhält man nach Elimination von R, und 4;: 
Xi—(k +1) Do Lkr takt!) 


2 2 k 5 en 
= Era o} ak layfı — (2 1% 1 eu A, . AR e RED au 1 [ax + ro Hr} 
(5,12). 


Beschränkt man sich beiderseits auf die niedrigste Potenz in r, (unter der Voraussetzung, daß 
1, „genügend‘‘ klein ist), so ergibt sich für 


k>1: XZurkHakapı. es SER She. (9,12; 
k=1: X1—-20, +3)" n4ıa— rad, X1— a)” r,a,a, (5,14). 
15 
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Fordert man in dem Bereich ® der Näherungen &, und der Lösung %, daß 


| s Mi, 

so ist: 

X sh MM Mh, Kk>l 
oder 

x —— | <a — ze M# Myyı . (5,15). 
Speziell ist also für : | 
BD: #—2|:<|o, —=! MM EREESMR ER. + U, TON 
k=3: 2— 2) le El MI am6,17); 
undsunk—l]: we E 

e-:1-Ml so -2P MM... 0.0.0... (5,18). 


Bei einmaliger Anwendung eines JV7ı (4,3) ist der Fehler 
—%|s |» — #2 M;. 
Vergleicht man diese Abschätzung mit (5,18), so sieht man, daß die Größenordnung in beiden 


Fällen gleich ist, d. h. die Verbesserung eines JV, aus drei Werten entspricht einer einmaligen 
Anwendung eines JV7r- 


Beispiel: 
Aus der quadratischen Gleichung & +ax +b=0 kann man das JV folgern: 
—b 
== ee Rd 
%n+1 = &, ee a = 
Mit a= — 5, b= 3ergibt sich die Iterationsfolge: 
U 
% —=U0,6, 


%, — 0,681 818, 
% = 0,694 737. 
Aus %,, %, %, folgt der verbesserte Wert & = 0,697 159. Fehlerabschätzung: M, = 0,162, 
M, = 0,038, &, — % = 0,097, also 
x —&| <7-10°. 
Der wahre Wert ist x = 0,697 224. 


Eine Wiederholung des Verfahrens aus 3 verbesserten Werten %,), &,, &, liefert einen 
neuen Wert mit einer Genauigkeit von 8 Ziffern. 


1,6. Verwendung iterierter Funktionen 
Bei einem JV @u+ı = f(z,) kann man das (p + 1)-te Glied der Iterationsfolge: 


Ins Entis +5 Into = Lpn+1s Intp+ls >» +5 Entop = Plnt2s::: =. . (6,1) 
aus dem ersten Glied der Folge durch das JV 
Pan = I ln) 0 ae 
bestimmen, wobei ?f (x) die p-fach Iterierte der Funktion f (x) ist. Es gilt also z. B. für 
=2: = = HN) = I: ern 
Ps: RE DE 


Von besonderem Interesse ist die Ordnung des JV (6,2). Es gilt der Satz: 


Für ein Vz &n4ı = f (@,) ist das JV mit der p-fach iterierten Funktion ?f (x) 


P&n+ı = ef Ca) 
von der Ordnung A?. 


m Dieses JV konvergiert, wenn man mit x = 0 beginnt, stets zu der Wurzel mit dem kleinsten Absolutbe- 
trag h - 


11) Pf (x) lies: p-te f(e). 
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Beweis alFurpap=2,.k=123: 


7a) =-/SU@)=F@; 
(=), 
Fe) = AS), 
f” (@) re DE) +3" Ne (@) HF" w 
an N are U later: 
TR... +0’ RR Hs PER) =. OFT TH. 
Re OSLO AL FR) =. 4120) +. 
Für eine ER: % ist also für f(®) = & unter der BAER: 
aleines Al, ee eye 
B) eines IH: FO)=0 "RÜLO F"@)+0,...; 
=") FO; 
ee N EA er 
FR)... yERM)-0, FIR) LO. 
bh) Rür psy = 1, 2: 
Ya =) =) = FW 


Ta =-fFWf@f(@), 
Ta) =") +- 


und daher: 


in a 2 ne RE ir af) +-. 
Ina) 2 (2 

le ei. Senf la) > Euiaetuunuek 
la) =... ART ee et 
Ya)... +8 WFT“ nr ° Sr 
ER) =... +48’ Wfl)" “ 


Für » = ®gilt für ein JV,, ‚fi Bde u I 
En IE Li w.z.b. w. 

Besonders wichtig ist, daß bei einem JV, auch das JV für die p-fach iterierte Funktion nur 
von 1. Ordnung ist. Der Vorteil der JV höherer Ordnung kommt auch dadurch wieder 
klar zum Ausdruck. 

Aus einem JVı +1 = f (&,) sollen unter Benutzung der iterierten Funktionen f(x) und 
5/ (x) Verfahren höherer Ordnung hergeleitet werden. 

a) Geht man aus von der Taylorreihe 


Ta) =) + af a) + a a" a) +: 


und führt wie in (5,3) das Differenzenschema der Iterationsfolge ein: 


% 
0 = Bw 


T 
i Gebt 


6@) a: Ö, N ö, | A (6,6), 
ER N | 
%g 
so ist für 
)=0,=f(%); = F(o) = (a): 
; 1 n 
6, = df (&) + Dr Ola) re ee (6,7); 


2) m= = (a), = Flo) = FR): 
ö2 +61 = (61 + du) f’ (Ro) F (dı +) +.» - 


15* 
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Zieht man (6,7) davon ar ist: 
1 [2 | 
PERL N ICh LT CHE ECT LUC HE PErEEE LE (6,8). 


Aus (6,7) und (6,8) folgt: 
&) Bricht man (6,7) nach dem ersten Glied ab, so ist 
Il) =. 5 ENG Aare en ee ER (6,9). 
0 « 
Die Anwendung des Ne wtonschen Verfahrens auf ® (x) = f(x) — « führt bei Ersatz von 
f(x) nach (6,9) zu einem JV &4ı = F («,) mit 


6 


F (a) = % — zB; 
0 


s ke Ve) Pe (6,10). 


Dies ist aber die in (5,5) hergeleitete Formel für den verbesserten Wert & aus drei Werten 
einer Iterationsfolge. Aus dieser Herleitung folgt auch sofort, daß das JV von 2. Ordnung ist 
und daß das Ausgangs-JV nicht zu konvergieren braucht. 


ß) Das Abbrechen der Taylorreihen (6,7) und (6,8) nach dem zweiten Glied liefert: 
pen (ö1 +2 60) — 60 0a Ög ö2 — 65 j 
0, 61 (do + 91) 0,01 (do + 81) 


Mit diesem Wert für f’(x,) kann man für ® (a) =f(«) — x das Newtonsche Verfahren 
2. Ordnung (3,7) umschreiben in die folgende rechte Seite des JV: 


5 055, (0 +8) E 
leerer See 


und Fis)=2 


str Koi 


Das Verfahren enthält statt der Ableitung die vier Werte der Iterationsfolge %, %1, %3 %z- 
Unter Benutzung von obigem f’(x,) und f’’(x,) lautet das Ne wtonsche Verfahren 
3. Ordnung (3,11): 


6, 62 — 61 
Kata ya) Tee aa 


b) Mit dem Ansatz für ein JV m 1 =F (®,): 
Fa)=/la) +AKFR-)+AKAFD- FR) >: (6,13) 


mit den freien Funktionen A, (x) und A,(x) lassen sich JV zweiter und dritter Ordnung 
aufstellen. 
Bei einem JV,; müssen erfüllt sein: 


Fe: unds rle)=0. 
Die erste Bedingnug ist erfüllt, die zweite ergibt 
FA-FO+HAAFO-YHALAFO-F@)=0. 
Da ?f’(#) =f'”(@), ergibt sich, wenn man fordert, daß diese Bedingung identisch erfüllt ist: 
AD -N HAAS - NH) =0 (*). 
Daraus läßt sich folgern: 
BRIRLEN: 
Lei 


BE f (@) er 
und die rechte Seite des JV lautet: 


Fo)=%+ Mao) 7 ne 1 (&, — %) 7 m Ss (% — #)|, (IV) (6,14) 


mit der willkürlichen Funktion A (x). 
Für ein JV 7], muß außerdem noch F”’ (&) = O sein. Man erhält 
F-F+2ARSFRD)—- 1) +24@(f@)- 1) 
TARA OT HN=0. 


„ii 
Pr 
Er 


Es 
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Zieht man davon die nach differenzierte Identität (*) für &=% ab und fordert, daß diese 
Bedingung identisch erfüllt werde, so ergibt sich: | 

Ad) AR 3A HN)HI=0 (**). 
Aus den beiden Identitäten (*) und (**) folgt: 


Aa=-ITI, Ald)=ra 


wa 4 = 1 


und die rechte Seite des JV 


An FOR F()) Fe) — /(a) 
Fa) =f(e) + ar A) ae: JVım) - - - (615). 
Bezeichnet man noch mit x; = f’ (a,), so ist nach (6,15) 
6, — 27 (2 —21)6, 
(1 — 1)? ? 


F(&)=% + (PVP ae (6.10): 


Dies ist ein JV777,, das nur einen ersten Differentialquotienten und drei Funktionswerte enthält. 
Beispiel: 
Aus der algebraischen Gleichung ® (x) = x° — 10 = 0 erhält man als eine rechte Seite 


eines divergierenden JV;: f (x) er und damit die Iterationsfolge: 


| 0%, = #0,167974 
x, —= 2,267 754 6, = — 0,322 764 
%; = 1,944 810 04 == .0,699.095., 
ü; = 2,643 903 
Es ergibt sich ’ 
nach (5,5) & = 2,15 727 (JV;r) mit einem Fehler von 3,4- 10°? ; 
(6,11) z= 2,15261 (JV7,) —1,8-10°3 ; 
(6,12) ©2153 955 (JVın) 5210787 
(6,16) & = 2,154 257 (J V 1) 1,8.--10>% 


Die Erhöhung der Ordnung eines JV um eins bringt den Gewinn von 1 Dezimalstelle. Die 
Ersetzung der Differentialquotienten durch Funktionswerte iterierter Funktionen ist wichtig 
bei solchen Funktionen, deren Ableitung nur umständlich zu bilden sind (z. B. bei nur tabuliert 
vorliegenden Funktionen). 


1,7. Iterationsvyeriahren’höherer Ordnung durch Linearkom- 
bination von Iterationsverfahren niederer Ordnung 


Es mögen eine Anzahl von JV; für dieselbe Wurzel # einer Gleichung gegeben sein: 


“o = fi (a0), (el URS EN ERRITTZTN: 
Da die Verfahren genau von %k-ter Ordnung sein sollen, gilt: 
Ra)=E ee Seel.) 


Es soll nun ein neues JV von mindestens (k + 1)-ter Ordnung gebildet werden durch Linear- 
kombination der Verfahren (7,1) mit zunächst freien ‚‚Gewichtsfunktionen‘‘ A;(x), dessen 
rechte Seite sei: 


p 
F(«) =2.A@) fa) ne. (7,3). 
Soll das neue Verfahren nun von (k + p — 1)-ter Ordnung sein, muß gelten: 
F@)=i, F@)=F'@)=..--Fetrm@)=-0, Fitrm@)z20 .. (74). 


Es git der Satz: 
Aus pJV k-ter Ordnung %4ı = fi(%), @=1,  - -, p), läßt sich, wenn die f; linear | (7,5) 


unabhängig sind, stets ein Verfahren (k + p — 1)-ter Ordnung herleiten. 


Be weis: Die »-te Ableitung von (7,3) lautet: 


Fofg) = en 0) AP (2) f® @): Ver (6). 


en) 


12)Hartree [4] weist auf diese Möglichkeit bin, glaubt aber mit Konstanten als Koeffizienten aus- 
zukommen. 
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Nach (7,2) und (7,4) muß gelten: 5 
FüÜ)=534@=8. 
i=1 


Fordert man die identische Erfüllung der letzten Gleichung, so gilt: 


53 Ala)=1... (77) und daher auch 2 A en (7,7%). 
Für die v-te Ableitung an der Stelle « = #gilt: 
ee FW(&) = % Z4P@) —Er 
Diese Bedingungen sind wegen (7,7*) erfüllt. 
b)vzk: 


a2, )araspa +, ara +- 


+40) =0. 
o)v=k: Aus 
Fo@) = 2 A@)R@)=0 
folgt die weitere Identität: =; 
SA) Tree (7,8,). 
P$»=kH+1: Aus = 
FED) = ZH) ASP) + A)NETE) = 0 


und durch Differenzieren von (7,8,) erhält man 


? 
> A;(e) ee) (Plz. 20 0 ae (7,3: +:). 
—ı 
In gleicher Weise kann man weitergehen und erhält: 
? 
SAFE E00) Se (7,8). 
i=1 


Um die Funktionen A;(x) zu bestimmen, muß man bis zur (k + p — 1)-ten Ableitung gehen. 
Die Identitäten (7,7) und (7,8) bilden ein lineares Gleichungssystem für die A;(z), das stets 
dann auch lösbar ist, wenn die f; linear unabhängig sind. Das neue JV ist daher genau von 
der Ordnung (kk +p — 1), w. z.b. w. 

Einfache Spezialfälle sind: 

(1) Aus zwei JV, ergibt sich ein JVr+ı (p = 2) mit der rechten Seite: 


(k) (k) 
B Oh | (U re ae 


(2) Aus drei JV, ergibt sich ein JV,,r (p = 3) mit der rechten Seite: 


ar: ee 
fe fs 8/1 2 : nd 
en + J. + er Ri mit D= 2 CENT) 


Speziell kann man diese Linearkombination von JV in der Weise vornehmen, daß man einer 
gegebenen Gleichung verschiedene Formen gibt und auf jede dieser Gleichungen z.B, das 
Newtonsche Verfahren anwendet, also ansetzt für die rechte Seite des JV: 


Fax)=% + 2 Ai(a) re 2 


2(7.il) 


Durch Differenzieren erhält man die Bedingungen, aus denen sich die A;(z) bestimmen lassen: 


Ikl@+t1=0, hf) = Zug) = ea gleir 


.. 


P 
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Für ein JV rn muß man daher von zwei Verfahren und den ersten beiden Gleichungen von 
(7,12) ausgehen, man erhält: 


Bla), () — Di(e) B;(a) 

Da) ©) ©i(a) 5 (a) 

In gleicher Weise folgt aus drei Verfahren ein JVıv (s. Beispiele). 
Beispiele: 


Es gibt verschiedene Wege, mehrere JV gleicher Ordnung für dieselbe Gleichung aufzu- 
stellen. 


a) Man kann von (3,8) ausgehen und mit verschiedenen Funktionen (x), z.B. pı(8) =1, 
?:(2) = x, ansetzen: 


F(«)=%& 


alhleriee rd) 


FE TOLION 
ercrng 


x: — a—=0 für die Quadratwurzel läßt sich folgern: 


’ == 


(SM 2,3,097, 


b) Aus der Gleichung (x) 


Tafelder behandelten Iterationsverfahren (Teil) 


Zahl der Rechenschritte 


Hinw. 
= des Iterationsverfahren > spez., wenn 2 
V allgemein ®Polyn. ext 
I I| 1 = f(an) F n 11 
2 I | ?anyı = "f(&n), wenn n4ı = f (in) ein JVı 3 
IV JVır 2F zn 6,3 
IX JVıu 
3 I I = an (£n); wenn In = if (£n) ein JVı 3 6,4 
LEN EA = vu |? 2 
4 Ei Pfnzı = f (in), wenn &nzı = [(@n) ein on pF Ex 6,2 
21 ag = DI) SNowtönsches.Merfah 2 3,7 
Inzı = in — ln) ewtonsches Verfahren F+R+l N \ 
? (&n) D (an) 
6 DH | 2241 = F+F*+F}* +2 2(n+k)+1| 38 
ende) 3 
Se D (an) p) 1 39 
2 II Cn4ı In FL (er) = jr ® (2)) F + HF -- 2 N = » 
D (in) + A (&n) D? (Xu) n * 
= m — F** 14 |O9n-+k,+k,+3| 3,10 
8 Il In+1 Tn & (zn) = u (£) ) (&) F r F, Ir F + = Ar al 2 
3 D (&n) D’ (zn) D? (an) 3% 4 311 
9 III a) 28% (2,) F+Ff,+PR+7 N = D 
® (zn) D’ (2) 
10 III —= in — — ee ER: 5 3n +2 3.12 
Tre) — HD len) P (en) has 
>* D (an) 1. Schritt: wie 10; 
11 I | inı = in — Fr (20) Zu 2.u.folg. Schritte: 2n+l1l 3,13 
12 | IH | zu = m — ; 
7 Tr) + Dam) Dem) , Pen | rm +m+7 I3n+4 | 3,14 
D’%2) +(k— 4) Dan) D’(an)  P’ (en) 
Se 1) 1. Schritt: wie 5; 
13 I | Znıı = &n — a 2.u.folg. Schritte: +1] @+1 [316 
0 
37 1. Schritt: wie 9; 
D (x) ©’ (x,) itkag 
14 Il nı= m -— Fe NZ, (a) a) 2. u. folg. Schritte 2n-+3 3,17 


F+R+4 
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es Zahl der Rechbenschritte Hinw. 
Im des Iterationsverfabren u spez., wenn Eu, 
Erz allgemein B Polsn. ex 


leerer — 
In+ı = in [6 8’2(zn) BD’) + 2 ® (&n) 8’ (&n) 
15 IV — 3 82 (24)] D(xn)/[6 D’2(&) D’ (zn) F+F+R+F+1l 4n+5 3,15 
+2®'(x,) D’’ (zn) — 6 D’’?(an)] ®’ (x) 


Dan) DD’ (zn) D (zn) 


en, = mars 
ra a a re 2 F+AR+R+R+17 | 4n+ll [3,18 
r ®: (£n) 13 ®'’2 (zn) — 6 (® ) 
ER. 6 D’4 (x) D (21) F 
® D’ (x) D? 
Inıı = in — vn — ra 1. Schritt: wie 16; 
ki hr ( A or 2. u. folg. Schritte: 3n+7 [3,19 
% % Ar 
rue 13 Fed _ ® zo) F+R+F+10 
1871] 3 Hm 2752 3 6,10 
N er 27 +6 = 5,9 
63.6, (60 + 61) 
20 gi Nr = 3F +8 _ 6,11 
Uls=ut 505, 720)- Ga Fa) vi“ 
2 BL 
EL do d, — 62 3F +17 = 6,12 
52. (6, + 26,)— 8 (86,+ 2 7° 
er x 
Fa) LEE) Een NSRTFrre = 6,14 
er +A(20) 2 (X) a] (1 X) % (zo) = (X 2) 
Ole 
N ne 2FLM+6 = 6,16 
2 
Bar, x Zn) — Jı (® &n) 
24 II | Znı = tn) rl A nl en FO+FW+FÜLF@+LZ - 7,9 
er &) rk) # 
235 | +1] mu = Er a 7® Se = FÜ +FW+F@ LFD +3 _ 7,9 
u n 
Br y 1 : ’ 
26 111 | In4ı = le n Kt hc fa 
Er Zn +79 +F®) + 16 _ 7,10 
+ A nen ely,l, — An > 
=. PT RE 
= ®; (£n) ®, (&n) — o, (zn) ®, (2) 2 2 K 


und damit aus f, und f, nach (7,7) und (7,8) die rechte Seite eines JV 


2.02 
ao) z@a— 32)" 


VD» 
oder aus fj, f»» fs: 


8a? 
(15a — lan? +3xt)’ 
Mit «>= 10, ©, —=.3,1 findet man: 
für das JV,; % = 3,16 410 mit einem Fehler von 18. 10%; 
JVyrr &ı = 3,16 234 „ 7 = 06 
c) Für die n-te Wurzel kann man die folgenden 3 Gleichungen aufstellen: 
B,()= tr —a=(0, D,(s) = ar — a2, D,(z) = ar — a8, 


JVım): 


.-— 
. 
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Nach (7,12) gilt: En 
Ahıt+h=-—1, (n— 1)A, +@n—1)%,=0 
oder 


h=— Vai ESF - 


Die A, sind in diesem Falle sogar konstant. 
Mitn = 3, «a = 10 und x, = 2 lautet die rechte Seite des JV 
| 52° + 50 2° — 100 
959 
und ©, = 2, &, = 2,153, x, = 2,154 434 (auf 7 Ziffern genau!). 
d) Aus DB (x) = x? — 10 = Ofolgt nach (7,12) für ein JVıv: 
Fe any Beer Ma ER EN BI Ren TR 


a 


oder 

20 16 5 
=—— Ag — 9? As N 9 P) 
und für 2, =2 wird z, = 2,15 471 mit einem Fehler von —28 - 10°, 


1,8: Rechenaufwand für die behandelten Iterationsverfahren 
— Erläuterungen zur Tafelder behandelten Iterationsverfahren — 


Die vorhergehenden Abschnitte haben gezeigt, daß es sehr zahlreiche JV jeder Ordnung 
gibt. Um die Wahl des geeignetsten JV für ein zu behandelndes Problem zu erleichtern, soll 
ein Überblick über den Rechenaufwand gegeben werden. Als Maß des Rechenaufwandes wird 
die Anzahl der Rechenschritte für die einmalige Anwendung des Verfahrens gezählt, und 
zwar unter der Annahme, daß zur Zahlenrechnung eine normale Rechenmaschine mit Hand- 
oder elektrischem Antrieb ohne besondere Speicherwerke verwendet wird. Gezählt werden 
daher als je ein Rechenschritt nur jede Multiplikation, Division oder Wurzelziehung unter 
Weglassung der Additionen und Subtraktionen. Die Zahl der Rechenschritte bei der Be- 
rechnung eines Funktionswertes f(x) oder ®(x) wird mit F, die der k-ten Ableitung mit F be- 
zeichnet. Im speziellen Fall kann man also durch Einsetzen die genaue Zahl der Rechenschritte 
angeben. Ist D oder fein Polynom, so ist #—=n, F,=n-—.1 (s. vorletzte Spalte). Für will- 
kürliche Funktionen p(%), A(&) usw. sei die Zahl der Rechenschritte F*. 


Zusammenfassung 
Ohne Vollständigkeit anzustreben, ist eine große Reihe von JV angegeben und der Weg 
gezeigt worden, weitere aufzufinden. Welches Verfahren man im Einzelfall einschlagen soll, 
ist natürlich nicht allgemein zu sagen. Die Ordnung des Verfahrens soll nach Möglichkeit stets 
größer als eins sein, je höher desto besser. Aber der Rechenaufwand darf dabei nicht unange- 
messen wachsen. Die Ordnung des Verfahrens und die Anzahl der für eine geforderte Genauig- 
keit notwendigen Schritte richtig vor Beginn der Rechnung gegeneinander abzugleichen, ist 
nicht leicht, wenn nicht sogar in vielen Fällen zu kompliziert und zeitraubend. So sind z.B. 
zur Lösung einer Gleichung n-ten Grades nach dem Ne wto.nschen Verfahren bei einmaliger 
Anwendung 2n Rechenschritte, bei einem JVır (s. lfd. Nr. 10) 3n + 2, bei einem JVıy (s. lfd. 
Nr. 16) 4n + 11 Rechenschritte notwendig. Schon hier sieht man, daß je nach der Größe von 
n die Entscheidung verschieden ausfallen wird, ob man ein JV;; mehrfach oder ein IV 

bzw. JVıy nur einmal oder wenige Male anwenden wird. 
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Singularly Loaded Rectilinear Plates-I 
by B.D. Aggarwala in Kharagpur (India) 


Mit Hilfe einer Analogie aus der Wirbeltheorie wurde die Anwendung nicht konvergierender oder schwach 
konvergierender Reihen bei der Berechnung der Durchbiegung und der Scherkräfte geradlinig begrenzter 
Platten unter Binzellast ausgeschaltet. Geschlossene Ausdrücks für die Scherkräfte werden als elliptische 
Funktionen erhalten, und eine Anzahl neuer zahlenmäßiger Ergebnisse in Tabellen gegeben. 


With the help of vortex analogy, the use ef non-convergent or slowly convergent series has been eliminated 
in calculating ihe deflection and Ihe shearing forces of rectilinear plates under concentrated loads. Closed. 
expressions for shearing forces are oblained in terms of elliptic functions and a number of new numerical 
resulis have been labulated. 


A l’aide d’une analogie de la theorie de tourbillon l’application des series non convergenles el peu con- 
vergentes au calcul de la deflection et des forces de tondage de plaques rectilignes sous charge concentr£e est 
&liminee. Des expressions fermees pour les forces de tondage sont obtenues en termes de fonctions elliptiques, 
et un nombre de nouveaux resultats num£riques est donne dans des tableauz. i 


IIo auasıorum us BHxpeBoi TeOpuH ÖkLIO HCKIIOYEHO IIPHMEHCHHE HECXONAINUXCH HH CTA0O 
CXONAINAXCA PANOB IPH pacyere IPoru6a MH CPesbIBAImMHX CHI IPAMONIHHeÄHO OTPAHHYEHHEIX 
IAIACTUH TION COCPeNOTOYeHHOH HaTpysKofi. 3AaMKHYTEIe BEIFAKEHHA MIA CPE3bIBAIOIIUX CHI 
HOIYYaMWTCA B BUN SINIHNTHYeCKUX PYHKIUM; MaeTca TAÖNMYHAA CBOMKA HOBEIX YyCJIOBBIX 
Pe3yıbTaTOß. 


1. Introduetion 
The problem of the bending of a rectilinear plate under an isolated load can be easily 
tackled with the help of conformal representation. Timoshenko [1] has applied the 
method of images and solved the problem for an equilateral triangular and an isosceles right 
angled triangular plate, the load acting at the centre in the first case and at an arbitrary point 
in the second one. But the method becomes rather cumbersome as we pass on to more compli- 
cated boundaries. 


In the present paper we have developed a method, making use of conformal represen- 
tation, with the help of which it will be possible to tackle certain types of singular loadings 
even when the boundary is fairly complicated and apply it to some particular problems. 


In the case of a concentrated load acting at any point A of the plate, we identify p? w 
(where wis the deflection and 77 — 0°/0x2 + 92/dy?) with the stream function of a vortex at 
Ain a cylinder of the same cross-section as the boundary of the plate. This at once leads us to 
an expression for 7} w, the integration of which gives w. The expression for y} w in itself is 
very useful in as much as it gives closed expressions for the shearing forces Q, and Q, in terms 

y of elliptice functions. These forces may be calculated 
N OA=@, Ax to any desired degree of accuracy. 

The advantage of this method of deriving Q, 
and Q, may be better made out with the help of an 
example. To get the shear at the mid-point of OC (see 
fig. 1) in the case of a square, we have, using the 


Re expression for wgivenby Timoshenko[2] 

S sin TR, 
4P Aare. 

C VeTen N where o, is the side of the square. 


This derivation involves differentiation of a dou- 
bly infinite series. This differentiated series being non- 
uniformly convergent, it is not easy to justify the process of differentiation. Moreover, the 
series for Q, is slowly convergent and useless for practical purposes. 


Using an alternative series expression for wgivenby Timoshenko[3], we get 
P ET MET 
er erPA. 5% a 


The process of differentiation is not easy to justify because the series for Q,is not even 
convergent. However, to three decimal places Tg$ = .917 and Tg (m&) = 1.000 for m>2 


so that we may write 
%=(P/o)(.17’—(l—1+1-1+ ...)) 


Figure 1 
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Using the Cessaro limit for the series 
1—-1+1-—1+--- 
we get Q, = (P/w,) (‚417) while our method gives simply Q, = P/m and we have 
o,/n = .4173135 . 
Comparing our results (2.7) with those of Timoshenko[2], we find 


ME. NN. NITZ 
MUSINZTZSIN= SIE 


WzE ei 2 2 ON 
er Dee or ech Dee ee) 


The right hand side may be looked upon asthe Fourier series expansion of the left hand 
side. Apparently, it would be quite difficult to obtain this identity otherwise. 

We take the plate in the planez= x +iy. IE P=P (x, y) be the intensity of the load 
at the point (z, y), the equation of equilibrium gives 

— pP 
D 

where pi = 02/022 + 02/öy?, wis the deflection in the direction of the normal to the plate 
(taken positive upwards), and Dis the flexural rigidity of the plate. 

In the case of a load acting at a point A (z,, Y,), the above equation reduces to the form 


BE 0HLE EN LE) 


except at the point A. 
The boundary conditions, when the edges are simply supported are 


w=0 TI DE EEE ra (1.3), 
ER a a ET TA) 
and 
- 0 
[Partes-r. ee et a ee et ER . (1.5). 
6 


C being the boundary arc of the plate and v the outward drawn normal direction. 

From (1.2) it follows that 77 wis a harmonic function which has a singularity at A and 
by (1.4) it vanishes at the boundary. It follows that 77 w may be identified with the stream 
function due to a rectilinear vortex inside a cylinder of the same cross-section as the boundary 
of the plate. If ® be the corresponding potential function, u the strength of the vortex, (1.5) 
gives 

S- D(d/w) pt wds = [ — D(o8/os)ds—=—u, sotht u=— P/D.. (1.6). 
c c 


Now the Schwartz-Christoffel transformation that transforms the interior 
of the plate in z-plane, into Im& > 0 may be written as 


A| BE -Wa+, a a N 


A,, B,, being suitable constants defining the position of the plate in the z-plane, Z,’s the points 
on the real axis in the Z-plane that correspond to the angular points in the z-figure, and zn a,'s 
are the internal angles at these points. 


In the case of a triangle, we may alwaystake, = — 1,, =0,{, =1. 
Now if the point £, corresponds to A, (1.6) gives 
Na > 
Hip a une, ; G=&— tn 
or AH 
pi w See IE rekey) 


== 08 —— 
anD T—i—% 
Equation (1.8) is particularly useful for calculating shearing forces Q, and Q, and in sections 
2 to 6 it is applied to calculate these forces in terms of elliptic functions for the following 
boundaries: — 
. Rectangle. 
. An equilateral triangle. 
. An isosceles triangle containing an angle of 120 degrees. 
. A right angled isosceles triangle. 
. A right angled triangle containing an angle of 30 degrees. 


aPp$wm— 
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2. Rectangle x 
Let @,, i®,, be the sides of the rectangle where @,, is real and @; a purely imaginary 
quantity. We form the Weierstrassian elliptic function p (z) with 2w, and 2w, as 
real and imaginary periods respectively. (1.7) may now be written sd —$ (2, 95; 93) where 
= 60 N (mm +2nm)t, 9— 140 2 2mo,+2n)*. 
N, N, 
@1, @,, @, being connected by the usual relation, namely, ©; +@ + @; — De 
(1.8) now becomes 5 
Er p2z—h2% pz—h% 2.1) 
2 Be et u N ee er ee Pi . Pe SE VEw BEL. Sr 7 . . . 
MOSEEDETZE N, ( 
(For sake of brevity, we shall not denote values of g,, g, every time. Unless otherwise stated 
they will always be assumed to be the same as used in the corresponding transformation). 


Therefore 


> 0,0 BT Dr jugate . . . 22 
a ren FH, + complex conjugate (2.2) 


and similary 


0 wPp[ 9% p'2 : 
Q,=— 2 DT. Be E ern] + complex conjugate . . (2.3). 


2, being the point of application of the load. 
In the particular case when the load acts at the centre, we have fh = 4,9 m =®&, 


yo; = e 
po = — es) Ha) ta +et, say 


where cis real and positive. 
(2.2) now becomes 


pP 2 42 N 
8 4 | EZ = TEE, —; Fr ; + complex conjugate 
es 2 RR) 
a I A + complex conjugate 
and (2.3) gives 
— cP iR S 
Q, = ar R a er complex Con jugatee 7 ee ey). 
For points on OC where pz = hp, p’ z = — b’ z, we get 
Q, = 0 
| p'z N EEERTE EEN 
£ n |(b2z—- 8)?+ 0 
In the particular case of a square g, — 0, and taking 9, = 4 so that ce = 1, we get 
=0 
| 9’ (y, 4,0) | ee RL 
= |1+9°(y,4,0) 
Similarly for points on OA where p zand p’z are both real, we get 
u 6'(z) N ee a De 
ge Fan 
 |b@)-o)+e | 
which for a square gives 
Q: —0 
—p (2.4, 0 
en en. 
= Ab om ,, Say 


(2.9) are of the same form as (2.7). 
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The values of « determining the shear in (2.9) for some values of & are ‚given in Table I. 


Tablel 
% 0 1/4 0/3 9/2 
& 0 ‚32180 41269 & 
Further we have for a square, so that , = (0, , —4, 
9’ (w, — ®) 26x 1 RE 


— [2000.07 


1a, — 2) (ha—1)} 
[p «— 17 

This relation illustrates the physically evident fact that values 2 shears at points of 0A 
equidistant from O and A are equal. 


3. Equilateral triangle 
In this case (1.7) becomes 


2= A, SE (@— 1) dc+ 3, 
we putö = 1/f, h®—= 1 — f?and adjusting the constants, we have [4] 
kh=p(l), %=0, G:=4 
24 
= Dr 
We take the angular points P,Q, Ras 2: ;o,/V3, 0, (1+ «/V3) respectively (see fig. 2): 
Therefore (1.8) becomes 


i 2% 
pyiw = 7; io g z El 2 —+ complex conjugate 
920 
at si + complex conjugate (3.1) 
res ATi an P J g EEE el Ei a & 
z, being the point of application of the load. 
Consequently we have 
P Hr 2 „ 2 . 
Q, ee u -rcomplex conjugate”". .ı= 2 (8.2), 
pP . „ Lad 2 e 
= - | 5 - = u = m a - complex conjugate -. . . . (3.3). 
Be 0 


In the particular case when the load acts at the centre, so that p’ z2,— Y12, we get 


P [73 2 „ 2 k 
„= — — — —————— | + complex conjugate 
ee = ZN aha | 
en 2 
ER 6 ee R | ankeomplex.conjugate Uwe au... 2 wu. (3.4) 
and er 
, 2 
A & nn H r - — 1 “RcorunleR- COMIUBaler!. Am a (3.5). 
For points on Q P where $ (z) is real, we get 
tn 
0 .—3J3P $(y,0,—4) 
zE 7 a AA ee ee (3.6) 
31302 
= —————g, say. 


IT 
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The values of a for a few values of y are given in Table 2. 


Table Re 
0 @j4y3 @/2y 3 auly 3 
. 0 .03072 .12201 1/3 
Bysymmetr , shears at other points may be calculated. (LhR 
a ne ns on the median through the origin, where z— — e?r, r denoting the distance 


1 la 
from the origin, 1, e, e? being the cube roots of unity such that e = D we 1-+i y3) 


= PD? 


dr 
a ER Be A NE (3.7) 
=: 2 per — 4 
and 
—9P p®r 
ee a rue. se a re 3.8 
d = an = 1 2 
Therefore 
14 cr 2 
pP Q = Qecos 80° 4 Qysin 0 ER Bir 
(3.9) 
aa a 
7t 
Yy @ 
a 0 X 
0 = X R © ? 
a AP=2@,, LQ=@, 
Figure 2 Figure 3 


The values of afor a few values of r are givenin Table 3. 


Table 3 
r 0 o,/4 ©/3 @,/2 2 [3 3 0,/4 @ 
a 0 09244 1/6 45534 00 — .95603 — 1/3 
As another numerical example, we may consider the shear at the point z—= — eo, when the 
load acts at any point 2, = — e?r on the median through the origin. In this case we get 
6P 1 
ee 
2 en re Ta; SFT ER 
EL 28 en 
al) 
Q = Q, cos 60° + Q, sin 60°— (6 P/a)[p (— r)]!. . ... . (3.11) 
= 3 Palm, say. 
The values of «for a few values of r are given in Table 4. 
Table 4 
r 0 o,/4 9/2 30/4 [ON 


& 0 .02798 22708 .90537 oo 
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4. Right Angled Isosceles Triangle 
In this case (1.7) becomes 


DEAN EUR N ARE ER EEME , (4). 
Writingd = 1/f, 1 — f? = 1/h?, and adjusting the constants, we write 
h=pe) n=4 —=0 
= (Ap??2)/(p 2). 


The angular points in the z-plane are taken to be 2 @, w(1-+i), 0; @, denoting the real half 
“period of the function p(2). 


(1.8) now becomes 
up 6’ 2, Hp? 2 — p3/2 2,6’ 2 
ee ae ee RE 
OT DE per 
Q = P 3/2 pU2 7 pH’ z.p’ 2 — ph" ap, £ 3/2 HR 2 p 2 — pp, 
Ag pop ?z2z—pzptlz, 62,922 — 9 zp8lez, (4.2) 


+ complex conjugate 


4 complex conjugate .... . AH, 


and 
ano ann, phtin. aEreB en urn) 
Ar Gar — ph zptz 6 2, 9°? 2 — bh’ zp?/2z, (4.3). 
+ complex conjugate | 
In particular, when the load acts on the axis of symmetry, so that 92, = bz, bh’. = — zu 
we get for the values of the shear at z= o,, wherebz2 =b2 =1,pz=k'z—=(. 
EN TREE. 
44.P [6°22, . 
|| re ats 
De LI A ee ae ce (4.4). 
4 9’ (@ı — Y) 
— 4Pa 
= ———, say 
7 
The values of & for some values of yare given in Table 5. 
Table5 
Yy 0 @,/4 0/2 3 01/4 oO 
& 00 .681795 22754 .05390 .00000 
Equating (4.4) to the corresponding result already known [1], we may postulate that 
ee 
al ee zı Ar N IT 2, ©] 2 @, . . . (4.5). 
Dale, en 


5. Isosceles Triangle Containing an Angle of 120 degrees 
In this case (1.7) becomes 
2 = A, SER (ee — 1) dc+ B,, 

setting & = 1/f, R®? —= 1/(1 — f?), and adjusting the constants [4], we write, 

h=hbl@), %=0, 9=4 

se 2 De 12 2 

ix ==; p’ 2% 4 
.We take the angular points P, Q, Rin the z-plane as 2 @,, @,(1-+ i//3), O respectively, &, being 
the real half period (see fig. 4). 
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Proceeding as in section 4, we get expressions for Q, and Q, . The shears Q, and Q, at the 
point z— w, when the load acts on the axis of symmetryat 2, — @, +? yare now given by 


9, = 0 
le 
u TERN, 
—P [p (y, 0, — 4) — 2]? Bee - 
oe (% 0, 4) pl Deren ]ie 21 “ 


on using the relations p(w, +74) = (bl) +/py) — 1) and Hl y) = —p(y,0, —4), 
the relations being immediate corrolary from the more general results [5] 
22 &aı-e)a—e) 
pe@+o)=&a+ p(2)— e s 
and 
| BARARg, Ag) = AH (2 9 95). 
Q 


0 TRP=20, 2 
Figure 4 
The values of a for a few values of y are given in Table 6. 
Table 6 
Yy 0 @/4 y3 @/2 y3 3@,/4 y3 0 
& es) 5.26673 1.97776 0.6311544 0 


6. Right Angled Triangle Containing an Angle of 30 degrees 
In this case (1.7) becomes 


De 4, je (€ a 18 ee 1)-2/3 dc B, : 
A writing & = YYf, h® = (f—1)/(f +1), and 


adjusting the constants, we write [4] 


h=he) H=0 —A4 

ee 1—h?z 

1+pz' 

We take the angular points P, Q, R in the 


z-plane as 0, @,, @(1 +i/Y3), ®, denoting 
f; PA =, q * the real half period of &(z) (see fig. 5). 


Figure 5 Now (1.8) gives 
1—p°2, 

— P 1+h?z2 1+9% 

2. en 0 
vo Da 0 
1i+p!2 1+$% 

NER Ha— hir, 2 —p37, 
In D ehyopen zn, ee (6.1) 


+ complex conjugate 
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differentiating which we get the expressions for Q, and Q,, namely, 


Rr3926 3p?2% 

le ar + complex conjugate . . . . . (6.2) 
Q BIN 36h?zh’ z 36?z6' 
" Anl pn H2—Hn 


On the real axis where pz=bz, 4 2=h'2,2= z, we get 


| Heomplex conjugate Eee): 


Q. =0 
ade baha Pape oe Er or 
 2n | —pz, bin pz, 
and on the line OR wherepz =bz,0 z = — k’z 
Q, nd 
Q ER Upizbe2 er 29 Te (6.5). 
 2n|hz— hr, Hz—pz, 
To take a numerical example we consider the load applied at the point at which 9°, = — $, i.e. 
92, = te?. This point is given by 
2, Ye. 929485 311172 er er (60% 
In this case Ihe formulae become 
—6:P[pzkh'z . 
er erg A complexconjugate Tr 6.7 
Q, Ar I wen A + complex conjugate (6.7) 
6 P[p:zhz : 
— — DE EEE T . (6.8). 
Q, > ker? \ + complex conjugate | (6.8) 
For points on the real axis where bz =bz, bh zz = 2,2 —=# 
Qz —0 
3P Seh 
ER SB Su Ol a Ar 6.9 
& an $a+1l 2 
= —3Pa/n, say 
The values of «for a few values of zare givenin Table 7. 
arb es 7 
Rn 0 o,/4 01/2 3 01/4 © 
& 0 .00515 69448 6341249 0 
similarly for points on QR 
Q, —=(0 
and =#: 
) _—3Pi p2zh'z . . . . . . . . . . . . . (6.10) 
eu. n. zart 


where z= w, +iy. Putting Q, = 3P a/n the values of & for a few values of y are given in 
Table 8. 


I7a:b1e»8 
y 0 @,/2Y 3 @/Y3 
& 0 ‚68988 0 


7. Regular Polygons (Defleetion) 


In general the integration of (1.8) will give an expression for w at any point. But it has 
not been found easy to integrate this equation as it stands. However a further transformation 


27: 
— a DR 71 
rt? er 


converts the z-area into the inside of a circle of radius unity in i{-plane. We can now re-write 


(1.7) as Br 
m Y ” a.—1 
AR) dE+B,. 


t 


16 


Se : 2: Math. Mech. 
234 Aggarwala, Singulariy Loaded Reciilinear Plates-T _na’38Nr.0 du 1806 


For a regular polygon of n sides, (7.1) takes the form 
t } : 
2A A Tr ee (7-2) 
0) 


And if the vertex t— 1 corresponds to z = c,, ©, being the radius of the circumcircle of the 
polygon, (7.2) gives on integration 


a 
0 ur 
where 1 
;& j 22 2 
nel Be ee BE DE 10 DEREN BE Se ee 
ln ) Bir Norzem, ee) 
and u > 3% 
hy = 2 Yyarı 
ern 1 1 2 
— ee: n\—2/n = — on: —— 2 a Pe a er . 
fe me de 2 1 .) (7.5) 
f) 


We get k, = 1.7666, k, — 1.3110, %,— 1.1744, k, = 1.1177. As ntends to infinity, k„tends to 1 
and the plate becomes circular. 
(1.8) now takes the form 


ee RT 
In’the particular case when i, = 0, i.e. the load acts at the centre of the plate, we get 


ern 
Vi ""AnD g 


AP = 
2 un 
Aue 
— Pı[- -dz - 
1040|? II ar [nr 


up . u Ben = 
= japan [ae [i 


on substituting zin terms of tfrom (7.3) this gives 
— 168wD - - je 
_—y— —= 22logtt — 2 Ayn+ı BD Agn+ı 19? dt + complex conjugate | \ 


Integrating term by term and he out the series, we get 
„—16rwD RE 
ne z log t 21212 ya Oki run hal + the complex conjugate | i 


On substituting 1 = g e'® we get 


— 16rwD En en 2 CH < x pn | 
D = 27 logtt — Eu Sr ertetroos nano]. 
P=0 g= i 


To this result we can add a harmonic function and because w = 0 on the boundary, we get 
finally, 


—16rwD E — a a 
rn jr > —. oPrtant2cos(gn— pn) o 


giving 


in - Sp n+1 %qgn+1 ER 
ESS ereoson— me] 
where |p n — q.n| denotes the absolute value of pn — qn. 
The central deflection w,is given by 
—- bw, D 25 ws @nı 
Pre 2m (7 
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The values of the constants «, for a few values of n are given in Table 9. 
Table 9 

nn Eee 4 We ern 

An 1.0082931 1.0022639 1.0008178 1.0004696 

n 7 8 9 KE Sep 

on 1.0001718 1.0000917 1.0000527 | 1.0000320 . 
IEn= 3and «a be the altitude of the equilateral triangle we get from above 

DE IE RLEN N Be a a (2,8). 


which value is an improvement upon the corresponding one given by Timoshenko [1] 
namely, (Pa?/D) .00571. 
In the case of n = 4, if a be the side of the square we get 


os lPaiDy. 0 REIHE. - (79) 


which is again nearer to the actual value’of &, than the corresponding value Siven by! Timo- 
shenko [2], namely (P a2/D) .01158. | 


Case of a Circle 
As ntends to infinity, a, tends to unity, and we get, 
a=Pa/8nD 


which is no t the central deflection of a simply supported circular plate when the load acts at 
its centre. 
The reason is not far to seek. For the boundary conditions for ne supported edges are 


and j ER 
0w 1o0w o9w 
Mai fe ze|-0 8 


where M as usual denotes the bending moment, o,the Poisso.nsratio of the material, P, the 
radius of curvature and ds the element of length in the direction of the tangent. In the case of a 
rectilinear plate 0%w/ds? — 0 due to boundary condition (a) and e is infinite at every point 
except the angular points, where it is undefined. Thus (b) reduces to 0?w/0v? = 0 which com- 
bined with 9?w/0s? = 0 is equivalent to p7w=0, the form that we have used. But as the 
number of sides increases such that each side becomes of the order ds in the formula for the 
radius of curvature 0— ds/dı, o does not tend to infinity but takes the value of the radius of the 
circumscribed or inscribed circle. In this limiting case the boundary condition that we have 
used does not hold good. The condition to be used becomes 
uw co dw 


ee 
This method has previously been employed by Seth [6] to calculate deflection of a 
clamped regular polygonal plate under uniform loading. This difficulty does not arise there be- 
cause the boundary conditions in that case are 


ow ’ 
Se in a ee re 


which conditions remain unaffected as the number of sides of the plate increases. He, therefore, 
obtains the results for a circle from those of a regular polygon by making n tend to infinity. 


Appendix 
Following is the list of the values of h(z) used in the pr together with an outline of the 
method we have used to derive them. 
resetag > ig, =0 
16* 
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| Table A(l) 
z 4 1/3 1/2 201/3 3 01/4 
9(2) st 5.27451 2.41421 1.46789 1.24009 


From the relation [7] 


a: en, El Pe 
enn" Ver een ER ) 


where cos? 9 = (s— 1)/(s + 1) we get, putting s — 1, k(sin »/4) = ®, Y2. But [8] K = 1.854075 
so that log w, = .1176123. However, a better approximation to the true value of log w, Is 
found to be [9] 0.1176122. 

Now from the series 


1 => 28 220 S —A.1.2 
Pott a er Zz 


we get the values of h(w,/4), K(@,/3). 
The value of p (w,/2) is found from the relation 


p (01/2) = eı + [(eı — &) (ei — es)” 


by putting , = 1, e,—=0, e, = — 1,andisequalto 1 V2: 
To find p (3@,/4) and p (2 @,/3) we put z= — w,/4, — @,/3 in the relation 
9 He) a a RR LE (A. 1.3) 
which in our case is equivalent to 
(z+0,) = a ee ar .: (A. 1.4). 


[The values of o(@, + i@,/4), $ (wı + i @,/3) ete.may also be found if we observe thatp(o, + iy)= 
$. (ey) + 1][p.(# y) — 1]! while from the relation p(A 2, 47% 9, 4 9,5) =A""h (, 95, 93) we get by 
puting 7=1,,—=4,9,—=0; p(i,4,0)—= —p (y, 4,0) so that p (m + iv) =by— 1] 
( (+1), also $ (w, — Yy) = (p y-+1)/(p y— 1) from (A. 1.4). So that 
: 1 
pw +iy)= p (wo, — Yy) > 
In this way we find that 
Ta’bre AN(2) 
2 © +im/4k © +io/3 0 +1/2 ©, +2:0/3 © +31@,/4 
9 (z) .806395, 68125, ‚41421 .18959 .10718 
2.1 Case 2.20, = 080, 4 
In this case we have the relation [10] 


p 
5 1/4 
a en (A. 2.1) 
(4 s® — 4) Y(1 — A? sin? p) 
where ; ı , 
EN Hr BE k=sin 15° 
s+(Y3—1) 


The corresponding value of K may be had from the tables [9]. We have log k = .2036154. 
Further from the relation @, — k/31’4 obtained from (A. 2.1) by putting $— 1, we get 


logo, 1084333 1 (A. 2.2). 
So that, using the series 
15,228 
Rz) = a: au 7 HOT Re ER (A.-2.3) 
we get the values of p (w,/4), p (w,/3). } 
To get the value of $ (w,/2), we use the relation (A. 1.3) which in this case becomes 


et ie a (AD 
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and get 


bw/)=1-+Y3. 
To get b (2 w,/3) we use the relation [11] 
39 (ui 3) = — pl) + g/p? u 
putting uv=2mw,/3 and using the relation [12] p(2iw,/Y3) = 0, we get p (2 0/3) = 41%, 


The value of p (3 ©,/4) may be obtained from 
To find the values of p (6 ©,/23), plio, 


the relation (A. 2.4). 
/Y3), etc. we use the relation 


—,]. „4 „10 5: a 
He) = = + Ta EFE +0 (2%) BEAT N) 
The values in this case are as follows: — 
TableA 6) 
% $@) 2 9(@) 

@,/4 10.85171 io,/&Y3 225514 
@,/3 6.10724 i@,/3 3 = 183098 
@1/2 2.73205 i@,/2y3 =8.13572 
20/3 1.58740 2io,/3/Y3 — 4.57074 
30/4 1.30452 350,/4 Y3 — 3.60591 
@, 1.00000 io/V3 — 2.00000 
2i30,/V3 — 0.00000 


To conclude, I wish to express my deep 


sense of gratitude to prof. B.R. Seth M.A,, 


D. Sc., Head of the Mathematics Department, for his suggesting me the problem and kind 


help throughout the preparation of this paper 
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KLEINE MITTEILUNGEN 


Numerische Lösung des Randwertproblems 
der Potentialgleichung mit Hilfe von Loch- 
karten. 


1. Allgemeines. 

In Weiterführung einer Arbeit von M. Kormes!) 
wurde ein Iterationsverfahren zur Lösung des Rand- 
wertproblems der Potentialgleichung 


au | 2 
022 Oy2 


mit Hilfe von Lochkartenmaschinen entwickelt. Grund- 
lage der Methode ist das bekannte Differenzenverfahren, 
bei dem das ebene Kontinuum durch ein Netz der 
Maschenweite h und die Differentialgleichung durch 
eine endliche Differenzengleichung ersetzt wird. Es 
werden zwecks Erreichung einer größeren Genauigkeit 


a (A), 


1) Marc Kormes, Numerical Solution ofthe Boundary Value 
Problem for the Potential Equation by Means of Punched Cards. The 
Review of Scientific Instruments 14 Nr. 8, 1943, S. 248— 250. 


auf Anregung von A. Weigand ?) zur Ermittlung des 
Funktionswertes in einem Netzpunkt die Werte in 
acht Nachbarpunkten herangezogen. Dabei gehen die 
Funktionswerte mit unterschiedlichem Gewicht in die 
Rechnung ein. Der Wert der unbekannten Funktion u 
für irgendeinen inneren Netzpunkt P;; wird: 

(2). 


Im Gegensatz dazu hat Kormes nur vier Nachbar- 
punkte benutzt, deren Funktionswerte bei der Rech- 
nung gleiches Gewicht haben. Wenn man die Funktions- 
werte in den Randpunkten kennt und für die inneren 


u 020 tur t ur 
+ %, £41) 
+0,05 (u, 14-1 T %rı,s-ı 
+% 42417 %rı, 241) 


2) A.Weigand, Die angenäherte Berechnung ebener und ro- 
tationssymmetrischer Potentialfelder mit Hilfe des Differenzen ver- 
fahrens. Verlag Technik, Berlin 1953. 
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02 ASTRA UML WME KIT m 1 20 2122 23 24 25 26 2728 


I SITE 
EEE EEENSEIZELTTTT 
anna m a mn 85 6 0 > 
h: a SE ERKREBEKENKE MR TR NS m 110% ma 102 m m 1 1 Io 

I 203 205 BERLEET 1 201910 TEE wa 365 1% 133 033 |, 
04204 % 


(iz 
I 25 25 25 26 26 26 26 26 265 27 265 26 255 25 29 15 12 17 169 168 167 167 167 167 17 |... 


j; “ 0305 3% air MN HM RR HM 0 27525 23 1 IE 


06 30 

nr BL. Ei 6 365 78 38 375 355 325 245 

ns 2.) DELL 1 M 4505 4 45 Ä6 3 265 
| A: Ze 


ei Eee 
N AEREENEEENETZTTETTN 


Rn TR 45 3 


r | ENSBTTER | 
12 2 BE A. Bild 1. Die Koordinatenachsen vermitteln 


4 | 100 | die Bezeichnung der Netzpunkte. Die in 

2 a Bee den Netzpunkten angeschriebenen Zahlen 

9 9 sind die dort geltenden Funktionswerte 
erster Näherung 


75 100 


17 RER wo| 8 60 4 
75 


80 6040 
2 01577 pe Te Pe 
arten! 
Spalte 1— 2 Nummer der Näherung s 
3 Punktart Randpunkt:1, innerer Punkt: 2 


4— 7  Punktbezeichnung 
8—11 Funktionswert 
12 Art der Karte Leitkarte:1, Anteilkarte:2 
13—16 Anteilfür Punkt Nr.: 
17—18 Faktor 
19—22 Produkt 


Tafel 2 
Zulochende Angaben in den Spalten: Zu lochende Angaben in den Spalten: 
3 4—7 8—11 | 12 | 3 47 8-34] 12 
1 0000 0000 1 1 0228 0067 1 
1 0001 0000 1 : : : i 
1 0002 0000 1 
; 1 0700 0350 1 
. : ; n 2 0701 0350 1 
1 0028 0000 al 5 i 3 2 
1 0100 0050 1 
2 0101 0050 1 
2 0102 0050 1 & & i : 
e : ; 5 2 0717 0245 1 
5 1 0718 0200 1 
2 0127 0033 ıl ß x 
1 0128 0033 1 ; e 5 : 
1 0200 0100 1 1 1810 0500 1 
2 0201 0100 1 2 1811 0750 1 
2 2 s 3 1 1812 1000 1 
il 1814 1000 1 
2 1815 0800 1 
2 1816 0600 A 
> 1817 0400 1 
1 1818 0200 1 
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' Netzpunkte ein System von Werten als erste Näherung 


der Funktionswerte sinnvoll wählt und darauf Glei- 
chung (2) anwendet, so erhält man für jeden Netzpunkt 
eine zweite Näherung. Bei genügend häufiger Wieder- 
holung des Verfahrens konvergieren die Funktionswerte 
gegen die Lösung des Systems von Differenzen- 
gleichungen. Geht h gegen Null, so nähert sich die 
Lösung des Systems der Differenzengleichungen der 
Lösung der Differentialgleichung (1)3). 

Zur Mechanisierung des Verfahrens werden je ein 
System von Leitkarten und ein System von Anteil- 
karten hergestellt, in die Punktbezeichnungen, Funk- 
tionswerte und Faktoren eingelocht werden, und eine 
Arbeitsanweisung durchgeführt, auf Grund derer die 
Karten so durch die einzelnen Maschinen geschickt 
werden, daß die Funktionswerte einer Näherung auto- 
matisch aus denen der nächst niedrigeren erhalten 
werden. 

Das Verfahren von Kor mes macht die Herstellung 
je eines Arbeitsmanuskriptes für die Leitkarten und die 
Anteilkarten erster Näherung erforderlich und damit 
auch das Ablochen und Prüfen dieser beiden Karten- 
sätze von Hand. Dabei ist zu beachten, daß das An- 
teilkartensystem bei Verwendung von vier Nachbar- 
punkten etwa viermal so umfangreich ist wie das Leit- 
kartensystem und bei Heranziehung von acht Nachbar- 
punkten entsprechend achtmal so groß, wenn man das 
Verfahren von Kormes sinnvoll auf diesen Fall 
überträgt. 

Das hier vorgeschlagene Verfahren gestattet eine 
mechanische Anteilkartenherstellung direkt aus den 
Leitkarten. Damit wird die Anfertigung des Arbeits- 
manuskriptes für die Anteilkarten und seine Ablochung 
und Prüfung von Hand erspart, bzw. durch Ausführung 
einer die im folgenden Abschnitt beschriebenen 14 Pro- 
grammpunkte (3) bis (16) umfassenden Arbeitsan- 
weisung ersetzt. 

Die Hauptvorteile der hier beschriebenen Methode 
sind die Schnelligkeit, mit der die Ergebnisse bestimmt 
werden, und der vollautomatische Ablauf der Arbeit. 


2. Einzelheiten über das Lochkarten- 
verfahren. 


Im folgenden wird zu den einzelnen Programm- 


punkten von Fall zu Fall auf ein Beispiel %) 5) Bezug 
genommen. Bild 1 zeigt das Gebiet, für das die Lösung 
der Potentialgleichung gesucht wird. Die vorgegebenen 
Funktionswerte in den Randpunkten und die Funk- 
tionswerte erster Näherung in deninneren Punkten sind 
eingetragen. Die Bezeichnung der Netzpunkte geht aus 
den Koordinatenangaben am Rande hervor. Jeder 
Punkt wird durch eine vierstellige Zahl gekennzeichnet. 
Die ersten beiden Ziffern geben die Nummer der Zeile, 
die letzten beiden die Nummer der Spalte. Den Loch- 
kartenentwurf zeigt Tafel 1, einen Auszug aus dem 
Arbeitsmanuskript Tafel 2. 


Die einzelnen Schritte des Rechenvorganges sind fol- 
gende: 


1. Die bekannten Randwerte und die angenommenen 
Werte für die inneren Punkte werden in den Spalten 
3—12 von Tafel 2 abgelocht. Es entsteht ein Leit- 
kartensatz. 

. Die unter (1) vorbereiteten Karten werden geprüft. 

. Die Leitkarten für die inneren Punkte werden 
heraussortiert. 

4. Zu den heraussortierten Karten werden je acht Teil- 
vervielfältigungen hergestellt als Grundstock für 
die Anteilkarten, derart, daß die Angaben der Leit- 
karten 
in den Spalten 3 4—7 


wm 


s) R.Courant, Z.f. angew. Math. Mech. 6 (1925), 322— 325. 

4) Die Aufgabe wurde von A.Weigan d gestellt im Anschluß 
an die oben zitierte Arbeit 8. 25. s. Fußnote 2, 8. 237. 

5) Für die sehr entgegenkommende Hilfe bei der praktischen Durch- 
führung der Arbeit danke ich Herrn Lachmun d vom Staatlichen 
Zentralamt für Statistik. 


11. 


12. 


13. 


14. 


15. 


16. 


17. 


22. 


23. 


24. 


25. 
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in den je acht Anteilkarten _ 
in den Spalten 
erscheinen. 


3 13—16 


. Die Leitkarten werden aussortiert und in den Spal- 


ten 1—2 mit Ol seriengelocht. 


. Die verbleibenden Anteilkarten werden in acht 


gleiche Gruppen aufgeteilt. _ 


. In der ersten Gruppe wird in den Spalten 4—7 die 


um 100 verminderte Zahl der Spalten 13—16 
eingelocht. 


. In der zweiten Gruppe wirdin den Spalten 4—7 die 


um 100 vermehrte Zahl der Spalten 13—16 
eingelocht. - 


. In der dritten Gruppe wird in den Spalten 4—7 die 


um 1 verminderte Zahl der ‘Spalten 13—16 
eingelocht. t 


. In der vierten Gruppe wird in den Spalten 4—7 die’ 


um l vermehrte Zahl der Spalten 13—16 
eingelocht. ; 
Die unter (7), (8), (9), (10) vorbereiteten Karten 
werden in Spalte 12 mit der Zahl 2 und in den Spal- 
ten 17—18 mit 20 seriengelocht. 
In der fünften Gruppe wird in den Spalten 4—7 
die um 10l verminderte Zahl der Spalten 
13—16 eingelocht. 

In der sechsten Gruppe wirdin den Spalten 4—7 die 
um 10l vermehrte Zahl der Spalten 13—16 
eingelocht. 

In der siebenten Gruppe wird in den Spalten 4—7 
die um 99 verminderte Zahl der Spalten 
13—16 eingelocht. 

In der achten Gruppe wird in den Spalten 4—7 die 
um 99 vermehrte Zahl der Spalten 13—16 
eingelocht. 

Die unter (12), (13), (14), (15) vorbereiteten Karten 
werdenin Spalte 12 mit 2undin den Spalten 17—18 
mit 05 seriengelocht. 

Die Anteilkarten werden nach den Spalten 4—7 
sortiert, um sie in die Reihenfolge der Bezeichnung 
der Netzpunkte zu bringen. 


. Die Leitkarten werden mit den Anteilkarten in der 


Reihenfolge der Spalten 4—7 gemischt und die An- 


gaben der Spalten 3 und 8—11 von den Leitkarten 


auf die Anteilkarten übertragen. 


. Die Leitkarten werden entfernt, indem man nach 


Spalte 12 sortiert. 


. Die verbleibenden Anteilkarten werden durch den 


Rechenlocher geschickt, um die Produkte der 
Funktionswerte (Spalten 8—11) mit dem Faktor 
(Spalten 17—18) zu erhalten. 


. Die Anteilkarten werden tabelliert, wobei nach den 


Spalten 13—16 geordaet und die zugehörigen An- 
gaben in den Spalten 19—22 summiert werden, 
so daß man den Gesamtwert aller zu einem Punkt 
13—16 gehörigen Anteile erhält. Für jeden Ge- 
samtwert stellt der Summenlocher eine Leitkarte 
zweiter Näherung her, in der Art, daß Angaben, 
die in den Anteilkartenin den Spalten: 


13—16 19—22 
in den Leitkartenin den Spalten: 
4—1 8—11 


erscheinen. In den Spalten 1—2 der Leitkarten 
wird 02, d.i. die Nummer der Näherung und in 
Spalte 3 die Zahl 2 eingelocht. 

Die Anteilkarten werden nach Spalte 3 sortiert. 
Die Karten für die Randpunkte, die in Spalte 3 die 
Nummer 1 tragen, werden weiter verwendet. 

Die Spalten 3—7, 12—18 für die inneren Punkte 
werden eine Anzahl von Malen gedoppelt, z. B. 
25 mal. 

Jeder der unter (23) entstandenen Sätzen wird inden 
Spalten 1—2 seriengelocht, beginnend mit 03, 
weiter 04,05... 

Nun werden aus den unter (21) entstandenen Leit- 
karten zweiter Näherung die Funktionswerte 
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(Spalten 4—7) in die Anteilkarten dritter Näherung 
übertragen und die Produkte (Spalten 8-11) x 
- (Spalten 17—18) gebildet. 

26. Die Anteilkarten dritter Näherung werden mit den 
unter (22) heraussortierten Karten für die Rand- 
punkte gemischt. 

27. Jetzt schließt sich ein Kreisprozeß an; und zwar 
werden nacheinander im Wechsel die unter Punkt 


BUCHBESPRECHUNGEN | ä 


Prof. Dr.-Ing. habil. Ulrich Graf und Dr. phil. Hans- 
Joachim Henning. Statistische Metho- 
den bei textilen Untersuchungen. 
X + 2788. m. 71 Abk. Berlin/Göttingen/Heidelberg 
1952. Springer-Verlag. Preis geb. 33,— DM. 

Das vorliegende Buch ist ein weiteres Zeichen für 
eine erfreuliche Tendenz, nämlich daß das Anwen- 
dungsgebiet der modernen statistischen Methoden 
nicht mehr auf das Feld der Forschung beschränkt 
bleibt, sondern auf die industrielle Produktion über- 
greift. Es ist auch kein Zufall, daß dieses Buch den 
textilen Untersuchungen gewidmet ist, sind doch 
gerade in Textilbetrieben schon längst laufend Messun- 
gen vorgenommen und notiert worden, die förmlich 
nach einer Auswertung drängen. 

Obwohl die Verfasser auf Beweisführungen weit- 
gehend verzichten, stellen sie an den Leser doch er- 
hebliche Anforderungen, denn ihm wird eine große 
Anzahl von Methoden vorgeführt, die zwingen, sich 
mit einer Menge neuer Begriffe vertraut zu machen. 
Das wird nur durch eine sehr geschickte Anordnung 
des Stoffes ermöglicht. Nach kurzen Einführungen 
werden die Methoden jeweils an ein oder mehreren 
Beispielen aus den verschiedensten Bereichen der 
textilen Untersuchung erläutert. Diese Beispiele bil- 
den den eigentlichen Kern des Buches, weil in ihnen 
die neuen Begriffe mit dem vertrauten Arbeitsbereich 
des Textilfachmannes in Verbindung gebracht werden. 
Jedes Beispiel bringt neue Gesichtspunkte zur Gel- 
tung und vermittelt wertvolle Ergänzungen vom sta- 
tistischen wie vom textilen Standpunkt. 

Der Inhalt des Buches ist durch die Stichworte: 
Grundbegriffe, Graphische Verfahren, Prüfung von 
Mittelwerten und Streuungen, Streuungsanalyse, Mu- 
tungsgrenzen, Korrelation und Kontrollkarten ge- 
kennzeichnet und wird durch je 10 Zahlentafeln und 
Nomogramme ergänzt. 

Referent, der Mitglied eines Ausschusses für stati- 
stische Auswertungsmethoden ist und daher die 
Schwierigkeiten bei der Einführung statistischer 
Verfahren in die Praxis kennt, möchte dieses aus- 
gezeichnete Buch jedem Textilfachmann wärmstens 
empfehlen, damit sich dieser von der Tragweite der 
dargestellten Methoden überzeugen und es als Leitfa- 
den für seine künftige Arbeit benutzen kann. 


Dresden. G Opitz 


James Jeans, Physik und Philosophie. 
3008. Zürich 1953. Rascher-Verlag. Preis geb. 
16,— DM. 

Verf. dieses außerordentlich anregend geschriebenen 
Buches sagt im Vorwort: „Der Zweck dieses Buches 
ist sehr einfach zu umschreiben; es ist die Erörterung 
— und bis zu einem gewissen Grade die Erforschung — 
jener Grenzgebiete zwischen Physik und Philosophie, 
die, lange Zeit als Ödland betrachtet, durch die 
jüngsten Entwicklungen der theoretischen Physik mit 
einemmal höchst wichtig und interessant geworden 
sind. Das neue Interesse erstreckt sich weit über die 
Fachprobleme der Physik und Philosophie bis zu 


2 
(21), (22), (25), (26) beschriebenen Operationen 
durchgeführt, wobei die einzulochende Nummer 


Anteilkarten unter (25) jeweils um eine Einheit 
fortschreiten. Das Verfahren wird solange fort- 
gesetzt, bis zwei aufeinander folgende Näherungen 
in den Spalten 8—11gleiche Werte zeigen. 

Berlin. Eleonore Schwarz. 


Tag 


Fragen, die das menschliche Leben sehr nahe angehen, 
wie Materialismus und Willensfreiheit. Deshalb, hoffe 
ich, wird das Buch viele interessieren, die weder Phy- 
siker noch Philosophen von Beruf sind, und zu diesem 
Ende habe ich die Erörterungen so einfach wie mög- 
lich gehalten, ...‘“ Der Stoff wird in sieben Ab- 
schnitten behandelt: I. Was ist Physik und was ist 
Philosophie ? II. Wie erkennen wir ? (Von Descartes 
zu Kant, Eddington). III. Die zwei Sprachen der 
Wissenschaft und der Philosophie. (Von Plato bis 
zur Gegenwart.) IV. Das mechanische Zeitalter zieht 
vorüber. (Von Newton zu Einstein.) V. Die neue 
Physik. (Planck, Rutherford, Bohr.) VI. Von der 
Erscheinung zur Wirklichkeit. (Bobr, Heisenberg, 
de Broglie, Schrödinger, Dirac.) VII. Einige Pro- 
bleme der Philosophie. Verf. kommt zu der Ansicht 
„die modernePhysik verhält sich nicht völligablehnend 
gegen einen objektiven Idealismus ähnlich dem von 
Hegel.“ und weiter „daß Determinismus und Freiheit, 
Materie und Materialismus im Lichte unserer neuen 
wissenschaftlichen Erkenntnis neu definiert werden 
müssen. Ist das einmal geschehen, dann muß der 
Materialist selber entscheiden, ob die einzige Art von 
Materialismus, welche die Wissenschaft noch zuläßt, 
passender Weise als Materialismus etikettiert werden 
kann“. 


Dresden. Willers. 


Proceedings of the Western Computer 
Conference. 231 S. m. zahlr. Abb. New York 1953. 
The Institute of Radio Engineers, Inc. Pr. geh. 3,50$. 

Das vorliegende Heft enthält den ausführlichen 
Text aller Vorträge, die auf der vom 4. bis 6. Februar 
53 in Los Angeles abgehaltenen Tagung über Ziffern- 
Rechenautomaten, Analogiegeräte usw. und ihre Ver- 
wendung gehalten wurden nebst den dazu gemachten 
Diskussionsbeiträgen. Auf dieser Tagung waren fünf 
Sektionen gebildet. In einer dieser Sektionen fand 
eine ausführliche Diskussion über die Entwicklung 
wie über Vorteile und Nachteile der Analogiegeräte 
und der Ziffernmaschinen statt, ferner über die 
Automatisierung des Rechnungsganges und der- 
artiger Vorrichtungen, wie sie z.B. zur Selbst- 
steuerung von Flugzeugen dienen. Eine andere 
Sektion beschäftigte sich mit der Anwendung von 
Rechenvorrichtungen im Geschäftsleben insbesondere 
im Versicherungs- und Steuerwesen. Eine dritte be- 
handelte die Verwendung in der Luftfahrt. In den 
beiden restlichen Sektionen wurde über sonstige An- 
wendungen dieser Maschinen und vor allem über 
technische Fragen gesprochen. So wurden neben 
vielem anderen über neue Speichervorrichtungen 
(Pulvari, Thorensen), über magnetische Reproduzier- 
und Druckvorrichtungen (Sims), über richtige Ver- 
wendung der Analogiegeräte (Cahn), über die Ver- 
wendung von Differentialanalysatoren zur Lösung 
partieller Differentialgleichungen (Howe u. Haneman), 
über einen neuen kleinen Differentialanalysator 
(Nordsieck) usw. berichtet. 


Dresden. Willers. 


Die besprochenen und angezeigten Bücher sind durch den Buchhandel zu beziehen, 


Herausgeber und Hauptschriftleiter: Prof. Dr. Dr.h.c. Fr. A. Willers, Dresden. Verlag: Akademie-Verlag GmbH., Berli 

Fernsprecher Sammelnummer: 200386. Postscheckkonto: Berlin 35021. Besiell und Vorleuse dieses "Heftes: ae 

schrift für angewandte Mathematik erscheint monatlich: Bezugspreis: vierteljährlich DM 15,—. Zuzüglich Bestellgeld. Einzelheft DM 6,— 

Abbestellungen können nur bis 4 Wochen vor Quartalsende anerkannt werden, andernfalls wird das folgende Quartalnoch geliefert. Veröffent- 

licht unter Lizenznummer 1207 des Amtesfür Literatur und Verlagswesen der Deutschen Demokratischen Republik. Satz und Druck : Druckerei 
„Ihomas Müntzer‘ Langensalza (57 275 4022) (1). Printed in Germany. 


Näherung in den Leitkarten unter (21) und in den. 


